2. Fabry-Perot cavity leaky-wave antennas (FPC-LWAs)

Transformation Optics in 2006 by Pendry [2], metasurfaces have represented a privileged platform for achieving a considerable control of electromagnetic waves propagation. Electromagnetic cloaking and metasurfing, i.e., controlling surface or guided waves through tunable metasurfaces [3], are one of the most known applications of metasurfaces. However, in this Chapter, we focus on the application of metasurfaces for the realization of reconfigurable leaky-

In the microwave range, there exists numerous realizations of metasurfaces (see, e.g., [4] and refs. therein), but at terahertz (THz) frequencies (nominally comprised between 300 GHz and 3 THz [5]), very few designs are available. Nowadays, THz technology is recognized as one of the most promising and challenging area of research for a twofold reason: (i) on the one hand, the wide and interdisciplinary character of THz applications, spanning from molecular spectroscopy and astrophysics, to high data rate communications and high-resolution imaging, passing through security screening and drug detection [6]; (ii) on the other hand, the increasing availability of efficient THz sensors and sources [7] that have recently contributed to

Nevertheless, the demand of efficient THz antennas is still high [6]. Indeed, even though various solutions have been proposed (see, e.g., [8–11]), efficient realizations often require high fabrication costs and complexity [11]. To better handle the efficiency vs. cost/complexity tradeoff, leaky-wave antennas [12] have been proposed and experimentally demonstrated as valid alternatives for the design of efficient, low-cost, THz antennas [13, 14]. Such prototypes [13, 14] present many advantages with respect to previous THz antenna solutions, but their lens-like structure does not allow for a planar design, an attractive feature for micromachined packaging of THz systems [11]. In this regard, more recently, a specific class of LWAs based on Fabry-Perot cavities (FPCs), i.e., FPC-LWAs, has been proposed as low-profile, fully-planar, directive, and efficient THz antennas [15]. Even more interestingly, tunable materials such as graphene [16] and liquid crystals [17] have recently been employed for the realization of FPC-

This Chapter is devoted to the analysis and design of novel THz FPC-LWAs. In Section 2, the general properties and the radiating features of FPC-LWAs are briefly reviewed. Conventional and more advanced designs are presented, with a specific focus on the technological constraints that typically affect the design of FPC-LWAs in the THz range. In Section 3, the design of THz FPC-LWA based on a homogenized metasurface is presented. The choice of a fishnetlike unit-cell is motivated by its remarkably low spatial dispersion and the capability of achieving a high reflectivity with a feasible variation of its structural parameters. The THz FPC-LWA considered in Section 3 produces a considerably high gain, but it does not allow for pattern reconfigurability. In Sections 4 and 5, two reconfigurable THz FPC-LWAs are therefore presented. In Section 4, the capabilities of two different THz FPC-LWAs based on graphene are shown and compared. The tunability of the graphene impedance through a DC voltage allows for a wide reconfigurability of the pattern. In Section 5, nematic liquid crystals are employed as tunable materials in a THz FPC-LWA to enhance the control of the beam angle through the

LWAs showing beam-steering capabilities at fixed frequency [18–22].

application of a low-driving voltage. Finally, conclusions are drawn in Section 6.

wave antennas (LWAs) [4].

90 Metamaterials and Metasurfaces

considerably narrow the so-called THz gap.

Fabry-Perot cavity leaky-wave antennas (FPC-LWAs) are partially open waveguiding structures, which support cylindrical leaky waves that radially propagate outward from the source [12, 23]. In this class of structures, radiation occurs through the excitation of the fundamental leaky modes supported by the structure, which are forward fast waves. Interestingly, when the excitation is a horizontal dipole (either electric or magnetic), the fundamental pair of TE, TM leaky modes is excited, and an FPC-LWA may produce a directive pencil beam at broadside or a conical beam with the cone axis along the vertical x-axis (see Figure 1) [12, 23, 24].

The architecture of an FPC-LWA relies on a grounded dielectric slab (GDS) covered with a partially reflecting screen (PRS) [25]. The gain enhancement phenomenon of these structures was originally explained through a ray-optics interpretation based on the Fabry-Perot concept [25] in the 1950s. However, the radiating features of these structures are more conveniently explained under the frame of the leaky-wave theory [12]. Indeed, due to the generality of the leaky-wave interpretation, the PRS can take various forms (e.g., a homogenized metasurface [24, 26] as in Figure 1(a), a denser dielectric superstrate [27, 28] as in Figure 1(b), and a distributed Bragg reflector [29] as in Figure 1(c)) provided that is suitably modeled with homogenized effective materials (when the PRS is bulky, e.g., a metamaterial) or surface impedances (when the PRS is planar, e.g., a metasurface). Regardless of the type of PRS, the leaky-wave approach yields to the same design process. This has allowed for the introduction of a variety of FPC-LWA designs, especially in the microwave range.

Figure 1. On the left, a typical FPC-LWA example fed by a horizontal magnetic dipole oriented along the y-axis. Depending on the operating point, such an antenna is able to produce either a pencil beam at broadside, or a conical beam. (a)–(f) Illustrative examples of Fabry-Perot cavity leaky-wave antennas (FPC-LWAs). (a)–(c) Conventional FPC-LWAs based on (a) a homogenized metasurface, (b) a high-permittivity dielectric cover layer (substrate-superstrate design), and (c) a distributed Bragg reflector (multiple layers dielectric structure). (d)–(f) Reconfigurable FPC-LWAs based on (a) a tunable graphene sheet, (b) a tunable graphene sheet in a substrate-superstrate configuration, and (c) a distributed Bragg reflector made of tunable nematic liquid crystals cells.

#### 2.1. Radiating properties of 2-D LWAs

The radiating features of an FPC-LWA qualitatively depend on the type of excitation that is considered [23]. In this Chapter, the FPC-LWAs under investigation are always assumed to be fed by horizontal magnetic dipoles (HMDs). This choice is motivated by (i) the constraints dictated by THz technology, as will appear more clearly in Section 2.6; and (ii) the interest in achieving broadside radiation, a feature that is never possible with vertical dipoles [23].

As it has been rigorously shown in [23], an HMD source launches a pair of leaky waves, one TM and one TE, both characterized by a (complex) leaky wavenumber kz ¼ β<sup>z</sup> � jα<sup>z</sup> 1 , where β<sup>z</sup> and α<sup>z</sup> are the phase and attenuation constants, respectively. Under certain conditions [30, 31] (usually met by a properly designed FPC-LWA [26]), the fundamental TE, TM pair of leaky modes is sufficient to describe radiation from such FPC-LWAs. In this case, the TM leaky wave determines the E-plane pattern, while the TE leaky wave determines the H-plane pattern [23]. For a broadside beam, these two leaky waves have very nearly the same phase and attenuation constants, and hence an omnidirectional pencil beam is created. Its beamwidths Δθ on the E- and H-planes are given by:

$$
\Delta\Theta\_{\rm E} = 2\sqrt{2}\hat{a}\_{z}^{\rm TM}, \quad \Delta\Theta\_{\rm H} = 2\sqrt{2}\hat{a}\_{z}^{\rm TE}, \tag{1}
$$

2.2. Design rules for FPC-LWAs

broadside or a conical beam).

is therefore given by ηrαrad=αz.

determined by the properties of the PRS.

2.3. Circuit model for FPC-LWAs

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi

<sup>0</sup>ε<sup>r</sup> � ð Þ <sup>m</sup>π=<sup>h</sup> <sup>2</sup>

its relative permittivity. Using Eq. (3) (under the assumption that β^

k 2

condition of broadside radiation (i.e., <sup>θ</sup><sup>0</sup> <sup>¼</sup> <sup>0</sup><sup>∘</sup>

q

<sup>z</sup> <sup>¼</sup> <sup>β</sup>TE <sup>z</sup> ¼

following design rule

<sup>β</sup><sup>z</sup> <sup>¼</sup> <sup>β</sup>TM

In all PRS-based FPC-LWAs, the PRS is used to create a leaky parallel plate waveguide (PPW) region, and the leaky waves are leaky (radiating) versions of the PPW guided modes that would be excited by the source in an ideal PPW, which results if the PRS is replaced by a perfect electric conductor (PEC) wall. This point of view allows for a simple design formula for the thickness of the dielectric layer in order to obtain a beam at a desired angle θ<sup>0</sup> (either a

For m-th order TM and TE PPW modes, the radial wavenumber of an ideal PPW would be

Typically, FPC-LWAs operate with the fundamental (i.e., m ¼ 1) leaky mode pair under the

It should be mentioned that the previous equations are found under the hypothesis that the FPC-LWA is of infinite transverse extent. This assumption holds as long as the leaky mode radiates most of its power as it approaches the lateral truncation. By defining the radiation efficiency η<sup>r</sup> as the ratio between the radiated power Prad and the initial power Pin, it results

This furnishes a design rule for the lateral truncation. Usually the antenna size is designed so that the FPC-LWA reaches at least a radiation efficiency η<sup>r</sup> ¼ 90%, and thus Eq. (5) reduces to L=λ<sup>0</sup> ¼ 0:366=α^z. However, this design rule assumes the structure to be lossless. When losses are present, the efficiency of the structure should be scaled of a factor α^rad=α^<sup>z</sup> [33], where α^rad is the leakage rate produced by the ideal lossless structure. The efficiency in a lossless structure

Finally, the optimal location of the source depends on its polarization. An HMD source maximizes the peak power density when it is placed on the ground plane [12]. However, the source location usually has little effect on the pattern shape, since this is dictated by the leakywave phase and attenuation constants. In particular, the phase constant is primarily determined by the thickness of the dielectric layer, while the attenuation constant is primarily

Regardless of the type of PRS, it is convenient to characterize its electromagnetic behavior in terms of a PRS impedance (admittance) ZPRS (YPRS ¼ 1=ZPRS). Under these assumptions,

θ0 � ��1=<sup>2</sup>

<sup>h</sup> <sup>¼</sup> <sup>m</sup>λ0=<sup>2</sup> <sup>ε</sup><sup>r</sup> � sin <sup>2</sup>

<sup>η</sup><sup>r</sup> <sup>≔</sup> <sup>P</sup>rad=Pin <sup>¼</sup> <sup>1</sup> � <sup>e</sup>�αzL (<sup>L</sup> being the transverse size of the cavity), and hence,

L=λ<sup>0</sup> ¼ �ln 1 � η<sup>r</sup>

, where h is the thickness of the dielectric substrate and ε<sup>r</sup>

Terahertz Leaky-Wave Antennas Based on Metasurfaces and Tunable Materials

); thus, the design rule reduces to <sup>h</sup> <sup>¼</sup> <sup>0</sup>:5λ0<sup>=</sup> ffiffiffiffi

� �=ð Þ <sup>2</sup>πα^<sup>z</sup> : (5)

<sup>z</sup> ≫ α^z), one finds the

εr p . 93

: (4)

http://dx.doi.org/10.5772/intechopen.78939

where the (� ̂ ) identifies normalization to the vacuum wavenumber <sup>k</sup><sup>0</sup> (i.e., ^kz <sup>¼</sup> kz=k<sup>0</sup> <sup>¼</sup> <sup>β</sup>^ <sup>z</sup> � jα^z), whereas the superscripts TE(TM) refers to the corresponding polarization.

For scanned beam β<sup>z</sup> ≫ αz, the beamwidths are given by [23]:

$$
\Delta\theta\_{\rm E} = 2\hat{\boldsymbol{\alpha}}\_{z}^{\rm TM} \sec \theta\_{0}^{\rm TM}, \quad \Delta\theta\_{\rm H} = 2\hat{\boldsymbol{\alpha}}\_{z}^{\rm TE} \sec \theta\_{0}^{\rm TE}.\tag{2}
$$

where θTE <sup>0</sup> (θTM <sup>0</sup> ) is the pointing angle (measured from the vertical x-axis) with respect to the H(E)-plane, whose expression reads [26]

$$
\sin \theta\_0 = \sqrt{\hat{\beta}\_z^2 - \hat{\alpha}\_z^2}.\tag{3}
$$

Equation (3) clearly defines the different radiating regimes of an FPC-LWA: when β<sup>z</sup> < αz, the antenna radiates at broadside, otherwise it radiates a scanned beam. At the limiting condition, i.e., when β<sup>z</sup> ¼ αz, also known as leaky cutoff condition or splitting point, the antenna radiates the maximum power density at broadside.

From Eqs. (1)–(3) it is clear that, as the wavenumbers of the two leaky waves begin to differ, the beamwidths and the pointing angles become different on the principal planes. This typically happens as the scan angle increases. Consequently, a desirable property for a PRS is to show almost the same behavior for both the TM and the TE polarization, especially when the angle of incidence changes. In fact, this would allow for having the equalization of the TM, TE leaky wavenumbers, which in turn produces nearly circular conical scanned beams [32].

<sup>1</sup> For the sake of simplicity, the z-axis is here chosen as the preferred axis of propagation.

#### 2.2. Design rules for FPC-LWAs

2.1. Radiating properties of 2-D LWAs

92 Metamaterials and Metasurfaces

E- and H-planes are given by:

̂

<sup>0</sup> (θTM

H(E)-plane, whose expression reads [26]

maximum power density at broadside.

where the (�

where θTE

1

The radiating features of an FPC-LWA qualitatively depend on the type of excitation that is considered [23]. In this Chapter, the FPC-LWAs under investigation are always assumed to be fed by horizontal magnetic dipoles (HMDs). This choice is motivated by (i) the constraints dictated by THz technology, as will appear more clearly in Section 2.6; and (ii) the interest in achieving broadside radiation, a feature that is never possible with vertical dipoles [23].

As it has been rigorously shown in [23], an HMD source launches a pair of leaky waves, one

and α<sup>z</sup> are the phase and attenuation constants, respectively. Under certain conditions [30, 31] (usually met by a properly designed FPC-LWA [26]), the fundamental TE, TM pair of leaky modes is sufficient to describe radiation from such FPC-LWAs. In this case, the TM leaky wave determines the E-plane pattern, while the TE leaky wave determines the H-plane pattern [23]. For a broadside beam, these two leaky waves have very nearly the same phase and attenuation constants, and hence an omnidirectional pencil beam is created. Its beamwidths Δθ on the

<sup>z</sup> , <sup>Δ</sup>θ<sup>H</sup> <sup>¼</sup> <sup>2</sup> ffiffiffi

<sup>0</sup> , <sup>Δ</sup>θ<sup>H</sup> <sup>¼</sup> <sup>2</sup>α^TE

<sup>0</sup> ) is the pointing angle (measured from the vertical x-axis) with respect to the

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi β^2 <sup>z</sup> � α^<sup>2</sup> z

q

Equation (3) clearly defines the different radiating regimes of an FPC-LWA: when β<sup>z</sup> < αz, the antenna radiates at broadside, otherwise it radiates a scanned beam. At the limiting condition, i.e., when β<sup>z</sup> ¼ αz, also known as leaky cutoff condition or splitting point, the antenna radiates the

From Eqs. (1)–(3) it is clear that, as the wavenumbers of the two leaky waves begin to differ, the beamwidths and the pointing angles become different on the principal planes. This typically happens as the scan angle increases. Consequently, a desirable property for a PRS is to show almost the same behavior for both the TM and the TE polarization, especially when the angle of incidence changes. In fact, this would allow for having the equalization of the TM, TE leaky

) identifies normalization to the vacuum wavenumber <sup>k</sup><sup>0</sup> (i.e., ^kz <sup>¼</sup> kz=k<sup>0</sup> <sup>¼</sup> <sup>β</sup>^

2 <sup>p</sup> <sup>α</sup>^TE

<sup>z</sup> sec θTE

1

<sup>z</sup> , (1)

<sup>0</sup> , (2)

: (3)

, where β<sup>z</sup>

<sup>z</sup> � jα^z),

TM and one TE, both characterized by a (complex) leaky wavenumber kz ¼ β<sup>z</sup> � jα<sup>z</sup>

<sup>Δ</sup>θ<sup>E</sup> <sup>¼</sup> <sup>2</sup> ffiffiffi

whereas the superscripts TE(TM) refers to the corresponding polarization.

For scanned beam β<sup>z</sup> ≫ αz, the beamwidths are given by [23]:

<sup>Δ</sup>θ<sup>E</sup> <sup>¼</sup> <sup>2</sup>α^TM

2 <sup>p</sup> <sup>α</sup>^TM

<sup>z</sup> sec θTM

sin θ<sup>0</sup> ¼

wavenumbers, which in turn produces nearly circular conical scanned beams [32].

For the sake of simplicity, the z-axis is here chosen as the preferred axis of propagation.

In all PRS-based FPC-LWAs, the PRS is used to create a leaky parallel plate waveguide (PPW) region, and the leaky waves are leaky (radiating) versions of the PPW guided modes that would be excited by the source in an ideal PPW, which results if the PRS is replaced by a perfect electric conductor (PEC) wall. This point of view allows for a simple design formula for the thickness of the dielectric layer in order to obtain a beam at a desired angle θ<sup>0</sup> (either a broadside or a conical beam).

For m-th order TM and TE PPW modes, the radial wavenumber of an ideal PPW would be <sup>β</sup><sup>z</sup> <sup>¼</sup> <sup>β</sup>TM <sup>z</sup> <sup>¼</sup> <sup>β</sup>TE <sup>z</sup> ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi k 2 <sup>0</sup>ε<sup>r</sup> � ð Þ <sup>m</sup>π=<sup>h</sup> <sup>2</sup> q , where h is the thickness of the dielectric substrate and ε<sup>r</sup> its relative permittivity. Using Eq. (3) (under the assumption that β^ <sup>z</sup> ≫ α^z), one finds the following design rule

$$h = m\lambda\_0/2\left(\varepsilon\_r - \sin^2\theta\_0\right)^{-1/2}.\tag{4}$$

Typically, FPC-LWAs operate with the fundamental (i.e., m ¼ 1) leaky mode pair under the condition of broadside radiation (i.e., <sup>θ</sup><sup>0</sup> <sup>¼</sup> <sup>0</sup><sup>∘</sup> ); thus, the design rule reduces to <sup>h</sup> <sup>¼</sup> <sup>0</sup>:5λ0<sup>=</sup> ffiffiffiffi εr p .

It should be mentioned that the previous equations are found under the hypothesis that the FPC-LWA is of infinite transverse extent. This assumption holds as long as the leaky mode radiates most of its power as it approaches the lateral truncation. By defining the radiation efficiency η<sup>r</sup> as the ratio between the radiated power Prad and the initial power Pin, it results <sup>η</sup><sup>r</sup> <sup>≔</sup> <sup>P</sup>rad=Pin <sup>¼</sup> <sup>1</sup> � <sup>e</sup>�αzL (<sup>L</sup> being the transverse size of the cavity), and hence,

$$L/\lambda\_0 = -\ln\left(1 - \eta\_r\right) / (2\pi\hat{\alpha}\_z). \tag{5}$$

This furnishes a design rule for the lateral truncation. Usually the antenna size is designed so that the FPC-LWA reaches at least a radiation efficiency η<sup>r</sup> ¼ 90%, and thus Eq. (5) reduces to L=λ<sup>0</sup> ¼ 0:366=α^z. However, this design rule assumes the structure to be lossless. When losses are present, the efficiency of the structure should be scaled of a factor α^rad=α^<sup>z</sup> [33], where α^rad is the leakage rate produced by the ideal lossless structure. The efficiency in a lossless structure is therefore given by ηrαrad=αz.

Finally, the optimal location of the source depends on its polarization. An HMD source maximizes the peak power density when it is placed on the ground plane [12]. However, the source location usually has little effect on the pattern shape, since this is dictated by the leakywave phase and attenuation constants. In particular, the phase constant is primarily determined by the thickness of the dielectric layer, while the attenuation constant is primarily determined by the properties of the PRS.

#### 2.3. Circuit model for FPC-LWAs

Regardless of the type of PRS, it is convenient to characterize its electromagnetic behavior in terms of a PRS impedance (admittance) ZPRS (YPRS ¼ 1=ZPRS). Under these assumptions,

Figure 2. (a) 2-D section of an FPC-LWA and (b) its equivalent circuit model.

FPC-LWAs are conveniently studied with an equivalent circuit model. Once ZPRS is known, the transverse resonance technique [12] can be applied to the equivalent circuit model (see Figure 2) to derive the relevant dispersion equations for the TE and TM modes

$$\left(Y\_0^p + Y\_{\text{PRS}}^p - jY\_1^p \cot\left(k\_0 h \hat{k}\_{x1}\right) = 0\right) \tag{6}$$

deal with materials which allow for describing the PRS impedance with a scalar value: more precisely, one for the TE polarization and one for the TM polarization. Indeed, we will always assume a negligible TE-TM cross coupling, thus reducing the dyadic surface impedance to a diagonal matrix. The previous dispersion equations will thus hold, provided that the correct

Terahertz Leaky-Wave Antennas Based on Metasurfaces and Tunable Materials

http://dx.doi.org/10.5772/intechopen.78939

95

As manifested from Eqs. (1)–(3), the beamwidths and the beam angles depend on the leaky wavenumber [12, 36, 37]. As the leaky wavenumber of an FPC-LWA depends on frequency (as a solution of the dispersion Eq. (6)), the radiating properties of a conventional FPC-LWA also depend on frequency. Indeed, the frequency scanning behavior is a distinctive feature common to all LWAs (not only FPC-type LWAs) [12]. All the structures outlined in this Section 2.4 are therefore referred to as conventional FPC-LWAs in the sense that they do not employ any

A homogenized metasurface (HMS) is characterized by a periodic lattice of "subresonant" elements, i.e., with period p ≪ λ0. Note that throughout the Chapter, we assume to deal with a patterned metallic screen, treating the metal as a PEC. This is still a good approximation at THz frequencies, provided that one considers thick enough (thicker than the skin depth [38])

The hypothesis p ≪ λ<sup>0</sup> is also known as the homogenization limit [34] and, as long as it is fulfilled, it allows for describing the electromagnetic behavior of the surface impedance through the generalized sheet transition conditions (GSTCs) [35]. Under this formalism, the tensorial surface impedance Zs is usually expressed as a diagonal matrix (no TE-TM coupling

ance Zs ¼ jXs (no ohmic losses are expected for HMSs made by PEC), Xs being the surface reactance. Furthermore, when dealing with HMS-based FPC-LWAs radiating at broadside,

Interestingly, as long as Xs ≪ η0, it has been shown [26] that the directivity at broadside D<sup>0</sup> is inversely proportional to the square of Xs, whereas the fractional bandwidth FBW (defined as the range of frequency, normalized to the operating frequency, for which the radiated power density at broadside decreases by less than 3 dB) is directly proportional to the square of Xs

, FBW <sup>¼</sup> <sup>2</sup> ffiffiffiffi

where tan δ is the loss tangent of the dielectric filling (which is nonnegligible in the THz range).

We note here that in the expression of D<sup>0</sup> published in [15] a factor 0:5=ε<sup>r</sup> was missing, and is now correctly included in

εr p π

X2 s η2 0

<sup>s</sup> are modeled with a purely imaginary imped-

<sup>s</sup> <sup>¼</sup> <sup>Z</sup>TE

<sup>s</sup> ¼ jXs.

þ tan δ, (8)

3

<sup>s</sup> , ZTM

PRS.

technique to avoid their intrinsic frequency scanning behavior.

good metals such as copper, aluminum, silver, and gold.

there is no distinction between the diagonal elements, and thus ZTM

<sup>D</sup><sup>0</sup> <sup>≃</sup> <sup>π</sup><sup>2</sup> 2 ffiffiffiffi ε3 r p

η2 0 X2 s

is assumed) whose diagonal elements ZTE

3

Eq. (8).

surface admittances are substituted to Y<sup>p</sup>

2.4. Conventional FPC-LWAs

2.4.1. Homogenized metasurfaces

where the equivalent admittances in air Y<sup>p</sup> <sup>0</sup> and inside the slab <sup>Y</sup><sup>p</sup> <sup>1</sup>, p∈ f g TE; TM for TE and TM modes have the following expressions:

$$Y\_0^{\rm TM} = \varepsilon\_\tau / \left(\hat{k}\_{\rm x0} \eta\_0\right), \qquad Y\_1^{\rm TM} = 1 / \left(\hat{k}\_{\rm x1} \eta\_0\right), \qquad Y\_0^{\rm TE} = \hat{k}\_{\rm x0} / \eta\_{0^\circ} \qquad Y\_1^{\rm TE} = \hat{k}\_{\rm x1} / \eta\_{0^\circ} \tag{7}$$

where <sup>η</sup><sup>0</sup> <sup>≃</sup>120π<sup>Ω</sup> is the vacuum impedance, whereas ^kx<sup>0</sup> <sup>¼</sup> ffiffiffiffiffiffiffiffiffiffiffiffiffi <sup>1</sup> � ^k<sup>2</sup> z q and ^kx<sup>1</sup> <sup>¼</sup> ffiffiffiffiffiffiffiffiffiffiffiffiffiffi <sup>ε</sup><sup>r</sup> � ^k<sup>2</sup> z q are the normalized vertical wavenumbers in air and in the dielectric, respectively. The zeros of Eq. (6) represent the wavenumbers of the eigenmodes of the structure in both the bound β^ <sup>z</sup> > 1 (i.e., surface waves) and the radiative β^ <sup>z</sup> <sup>&</sup>lt; 1 (i.e., fast leaky waves) regimes<sup>2</sup> .

In the case of a substrate-superstrate, or multiple-layer configurations (see Figures 1(a) and (b)), the layers can be treated as transmission line segments, and the resulting PRS impedance can be modeled as the input impedance seen looking upward from the substrate interface. Such an input impedance will take two different expressions depending on whether the TM or the TE polarization is considered. The previous dispersion equations will thus hold, provided that the correct input admittances are substituted to Y<sup>p</sup> PRS.

In the case of homogenized metasurfaces, the concept of surface impedance is conveniently introduced (see [26, 34]). While in general the surface impedance is described by a tensorial surface impedance Zs (see [35] for further details), throughout the Chapter, we will always

<sup>2</sup> It is worth recalling that surface waves are proper waves (improper surface waves exist from a mere mathematical point of view, but they are not physically meaningful) since they exponentially decay in the vertical direction α<sup>x</sup> > 0, whereas forward (β^ <sup>z</sup> > 0) leaky waves are improper waves, as they exponentially grow in the vertical direction α<sup>x</sup> < 0. This does not compromise the physical significance of leaky waves, since the exponentially-growing character of leaky waves is maintained only within a limited angular region, and thus they do not violate the Sommerfeld radiation condition at infinity (see [30, 31] for a rigorous discussion).

deal with materials which allow for describing the PRS impedance with a scalar value: more precisely, one for the TE polarization and one for the TM polarization. Indeed, we will always assume a negligible TE-TM cross coupling, thus reducing the dyadic surface impedance to a diagonal matrix. The previous dispersion equations will thus hold, provided that the correct surface admittances are substituted to Y<sup>p</sup> PRS.

#### 2.4. Conventional FPC-LWAs

FPC-LWAs are conveniently studied with an equivalent circuit model. Once ZPRS is known, the transverse resonance technique [12] can be applied to the equivalent circuit model

> <sup>1</sup> cot <sup>k</sup>0h^kx<sup>1</sup> � �

<sup>0</sup> and inside the slab <sup>Y</sup><sup>p</sup>

, YTE

<sup>0</sup> <sup>¼</sup> ^kx0=η0, YTE

q

<sup>z</sup> <sup>&</sup>lt; 1 (i.e., fast leaky waves) regimes<sup>2</sup>

PRS.

ffiffiffiffiffiffiffiffiffiffiffiffiffi <sup>1</sup> � ^k<sup>2</sup> z

¼ 0, (6)

and ^kx<sup>1</sup> <sup>¼</sup>

.

<sup>1</sup>, p∈ f g TE; TM for TE and TM

<sup>1</sup> <sup>¼</sup> ^kx1=η0, (7)

q

ffiffiffiffiffiffiffiffiffiffiffiffiffiffi <sup>ε</sup><sup>r</sup> � ^k<sup>2</sup> z

are

<sup>z</sup> > 1

(see Figure 2) to derive the relevant dispersion equations for the TE and TM modes

PRS � jY<sup>p</sup>

<sup>1</sup> <sup>¼</sup> <sup>1</sup><sup>=</sup> ^kx1η<sup>0</sup>

� �

the normalized vertical wavenumbers in air and in the dielectric, respectively. The zeros of Eq. (6) represent the wavenumbers of the eigenmodes of the structure in both the bound β^

In the case of a substrate-superstrate, or multiple-layer configurations (see Figures 1(a) and (b)), the layers can be treated as transmission line segments, and the resulting PRS impedance can be modeled as the input impedance seen looking upward from the substrate interface. Such an input impedance will take two different expressions depending on whether the TM or the TE polarization is considered. The previous dispersion equations will thus hold, provided

In the case of homogenized metasurfaces, the concept of surface impedance is conveniently introduced (see [26, 34]). While in general the surface impedance is described by a tensorial surface impedance Zs (see [35] for further details), throughout the Chapter, we will always

It is worth recalling that surface waves are proper waves (improper surface waves exist from a mere mathematical point of view, but they are not physically meaningful) since they exponentially decay in the vertical direction α<sup>x</sup> > 0, whereas

not compromise the physical significance of leaky waves, since the exponentially-growing character of leaky waves is maintained only within a limited angular region, and thus they do not violate the Sommerfeld radiation condition at

<sup>z</sup> > 0) leaky waves are improper waves, as they exponentially grow in the vertical direction α<sup>x</sup> < 0. This does

Yp <sup>0</sup> <sup>þ</sup> <sup>Y</sup><sup>p</sup>

Figure 2. (a) 2-D section of an FPC-LWA and (b) its equivalent circuit model.

, YTM

where <sup>η</sup><sup>0</sup> <sup>≃</sup>120π<sup>Ω</sup> is the vacuum impedance, whereas ^kx<sup>0</sup> <sup>¼</sup>

that the correct input admittances are substituted to Y<sup>p</sup>

where the equivalent admittances in air Y<sup>p</sup>

� �

modes have the following expressions:

(i.e., surface waves) and the radiative β^

infinity (see [30, 31] for a rigorous discussion).

<sup>0</sup> <sup>¼</sup> <sup>ε</sup>r<sup>=</sup> ^kx0η<sup>0</sup>

YTM

94 Metamaterials and Metasurfaces

2

forward (β^

As manifested from Eqs. (1)–(3), the beamwidths and the beam angles depend on the leaky wavenumber [12, 36, 37]. As the leaky wavenumber of an FPC-LWA depends on frequency (as a solution of the dispersion Eq. (6)), the radiating properties of a conventional FPC-LWA also depend on frequency. Indeed, the frequency scanning behavior is a distinctive feature common to all LWAs (not only FPC-type LWAs) [12]. All the structures outlined in this Section 2.4 are therefore referred to as conventional FPC-LWAs in the sense that they do not employ any technique to avoid their intrinsic frequency scanning behavior.

#### 2.4.1. Homogenized metasurfaces

A homogenized metasurface (HMS) is characterized by a periodic lattice of "subresonant" elements, i.e., with period p ≪ λ0. Note that throughout the Chapter, we assume to deal with a patterned metallic screen, treating the metal as a PEC. This is still a good approximation at THz frequencies, provided that one considers thick enough (thicker than the skin depth [38]) good metals such as copper, aluminum, silver, and gold.

The hypothesis p ≪ λ<sup>0</sup> is also known as the homogenization limit [34] and, as long as it is fulfilled, it allows for describing the electromagnetic behavior of the surface impedance through the generalized sheet transition conditions (GSTCs) [35]. Under this formalism, the tensorial surface impedance Zs is usually expressed as a diagonal matrix (no TE-TM coupling is assumed) whose diagonal elements ZTE <sup>s</sup> , ZTM <sup>s</sup> are modeled with a purely imaginary impedance Zs ¼ jXs (no ohmic losses are expected for HMSs made by PEC), Xs being the surface reactance. Furthermore, when dealing with HMS-based FPC-LWAs radiating at broadside, there is no distinction between the diagonal elements, and thus ZTM <sup>s</sup> <sup>¼</sup> <sup>Z</sup>TE <sup>s</sup> ¼ jXs.

Interestingly, as long as Xs ≪ η0, it has been shown [26] that the directivity at broadside D<sup>0</sup> is inversely proportional to the square of Xs, whereas the fractional bandwidth FBW (defined as the range of frequency, normalized to the operating frequency, for which the radiated power density at broadside decreases by less than 3 dB) is directly proportional to the square of Xs 3

$$D\_0 \simeq \frac{\pi^2}{2\sqrt{\varepsilon\_r^3}} \frac{\eta\_0^2}{X\_s^{2'}}, \quad \text{FBW} = \frac{2\sqrt{\varepsilon\_r}}{\pi} \frac{X\_s^2}{\eta\_0^2} + \tan \delta,\tag{8}$$

where tan δ is the loss tangent of the dielectric filling (which is nonnegligible in the THz range).

<sup>3</sup> We note here that in the expression of D<sup>0</sup> published in [15] a factor 0:5=ε<sup>r</sup> was missing, and is now correctly included in Eq. (8).

As manifested in Eq. (8), a tradeoff is established between FBW and D0. Consequently, lower values of Xs lead to highly directive FPC-LWAs working in a narrow bandwidth, whereas higher values of Xs lead to moderately-directive FPC-LWAs working in a considerable bandwidth. To give some numbers, directivities higher than 30 dB can be obtained with Xs < 20 Ω, assuming ε<sup>r</sup> ¼ 2:3.

2.4.2. Substrate-superstrate structure

PRS acts as a more reflective surface.

2.4.3. Multilayered dielectric structure

2.5. Reconfigurable FPC-LWAs

As shown in [27, 28], a very high step index (i.e., ffiffiffiffiffiffiffiffiffiffiffiffi

upward from the substrate interface) of a factor ð Þ <sup>ε</sup>r2=ε<sup>r</sup> <sup>n</sup>

In a substrate-superstrate configuration, the PRS is represented by a cover layer made of a denser dielectric material with respect to the substrate, as depicted in Figure 1(b). These kinds of structures, commonly known as substrate-superstrate (SS) LWAs, have been extensively studied by N. G. Alexópoulos and D. R. Jackson [27, 40]. In those works [27, 40], they showed that the PRS (represented by the superstrate) would act as an almost perfectly reflecting layer, when the thickness h<sup>2</sup> of the superstrate (the cover layer) is chosen so that it is an odd multiple of one-quarter of the wavelength in the vertical direction. This corresponds to the condition

<sup>h</sup><sup>2</sup> <sup>¼</sup> ð Þ <sup>2</sup><sup>m</sup> � <sup>1</sup> <sup>λ</sup>=<sup>4</sup> <sup>ε</sup>r<sup>2</sup> � sin <sup>2</sup>

It is clear from Eqs. (4) and (9) that the design rules for optimizing radiation at broadside (θ<sup>0</sup> ¼ 0) with the fundamental (m ¼ 1) TE, TM leaky modes suggest a thickness of a half- and of a quarter-wavelength in their respective media for the substrate and the superstrate, respectively. Another distinctive feature of SS-LWAs is that the directivity increases as the permittivity of the superstrate layer increases with respect to that of the substrate layer since the superstrate

high directivities with a substrate-superstrate configurations. As high-permittivity materials usually exhibit also nonnegligible loss tangent, the radiation efficiency of a substrate-superstrate FPC-LWA is often limited by the relatively high dielectric losses of the dielectric layers. In this regard, a different way to increase directivity was proposed in [41]. There, the single superstrate is replaced by a periodic array of such superstrates, as shown in Figure 1(c). In this structure, the PRS consists of a stack of multiple alternating layers of high permittivity and low permittivity, whose thicknesses are chosen according to Eq. (9). As shown in [41], the directivity increases geometrically with the number of superstrate layers, and thus very directive beams may be obtained using modest values of superstrate permittivity, provided that several superstrate layers are used. From a circuital point of view, the gain enhancement phenomenon of this class of structures is readily explained by noticing that the 2n � 1 quarter-wavelength sections reduce the free-space impedance (seen as a load impedance of the transmission line segment looking

As a last remark, it is worth noting that the multistack of alternating layers can be interpreted as a distributed Bragg reflector (DBR) meaning that it is a periodic electromagnetic bandgap (EBG) structure operating in a stopband. From this point of view, the substrate acts as a defect in the periodic EBG structure [42]. More interestingly, a leaky-wave explanation of the phenomenon is provided in [29].

In most of microwave applications, the frequency scanning behavior typical of conventional FPC-LWAs is an undesirable effect, since the antenna radiating performance should not change within a small fractional bandwidth centered around the operating frequency. Conversely, it

θ0 � ��1=<sup>2</sup>

Terahertz Leaky-Wave Antennas Based on Metasurfaces and Tunable Materials

εr2=ε<sup>r</sup>

.

: (9)

http://dx.doi.org/10.5772/intechopen.78939

97

p ) is required for achieving relatively

The previous relations hold for broadside radiation. When the antenna radiates a scanned beam, the TM and the TE surface impedances of an HMS generally differ each other. Indeed, the dependence of the surface impedance on the angle of incidence (spatial dispersion) and frequency (temporal dispersion) should be taken into account in a thorough characterization of an HMS. A canonical example is that of patch array and its complementary structure (according to Babinet's principle), the strip grating (see Figure 3(a)). Indeed, the grid impedance of an array of patches (a strip grating) for TE(TM) polarization exhibits a spatially dispersive character because it depends on the angle of incidence θ. The formula for the HMS of a strip grating and a patch array of infinite extent are given in [34], and are here applied for obtaining results in Figures 3(b) and (c), assuming a lattice of period p ¼ λ=10 with w ¼ p=10 and g ¼ p=10 for the strips and the patches, respectively. As is manifest from Figures 3(b) and (c) the presented topologies (strip and patch) allow for synthetizing low impedances only for unfeasible values of their structural parameters (i.e., g and w) in the THz range. For example, at 1 THz and for <sup>θ</sup> <sup>¼</sup> <sup>0</sup><sup>∘</sup> , an impedance of 20 Ω is reached with a strip grating with strips narrower than 3 μm. As it will appear clearly in Section 2.6, this value could present some issues during the photolithography fabrication process [15]. Moreover, the high variation of the reactance with respect to the angle of incidence would produce very different radiating features as the beam is scanned from broadside to endfire.

However, several HMS geometries with limited spatial dispersion for both TE and TM polarizations have been proposed in the literature [35, 39], even for the development of highly directive FPC-LWAs [15]. Among them, the so-called fishnet geometry (see Figure 3(a)), i.e., a superposition between the patch and the strip layouts, will be investigated in Section 3. We note that the synthesis of HMSs with limited spatial dispersion is a key aspect for the development of reconfigurable LWAs. In fact, it allows for designing and optimizing the radiating properties of a reconfigurable FPC-LWA at broadside (i.e., <sup>θ</sup> <sup>¼</sup> <sup>0</sup><sup>∘</sup> ), without affecting the impedance properties as the beam is steered.

Figure 3. (a) Unit cells of three different type of HMSs: the strip grating, the array of patches, and the fishnet. (b) and (c) Surface reactance Xs vs. frequency in the range 0:5 < f < 1:5 [THz] for angles of incidence in the range 0<sup>∘</sup> ≤ θ ≤ 80<sup>∘</sup> with angular steps of 10<sup>∘</sup> : (b) strip grating for TM polarization and (c) array of patches for TE polarization (parameters in the text).

#### 2.4.2. Substrate-superstrate structure

As manifested in Eq. (8), a tradeoff is established between FBW and D0. Consequently, lower values of Xs lead to highly directive FPC-LWAs working in a narrow bandwidth, whereas higher values of Xs lead to moderately-directive FPC-LWAs working in a considerable bandwidth. To give some numbers, directivities higher than 30 dB can be obtained with Xs < 20 Ω,

The previous relations hold for broadside radiation. When the antenna radiates a scanned beam, the TM and the TE surface impedances of an HMS generally differ each other. Indeed, the dependence of the surface impedance on the angle of incidence (spatial dispersion) and frequency (temporal dispersion) should be taken into account in a thorough characterization of an HMS. A canonical example is that of patch array and its complementary structure (according to Babinet's principle), the strip grating (see Figure 3(a)). Indeed, the grid impedance of an array of patches (a strip grating) for TE(TM) polarization exhibits a spatially dispersive character because it depends on the angle of incidence θ. The formula for the HMS of a strip grating and a patch array of infinite extent are given in [34], and are here applied for obtaining results in Figures 3(b) and (c), assuming a lattice of period p ¼ λ=10 with w ¼ p=10 and g ¼ p=10 for the strips and the patches, respectively. As is manifest from Figures 3(b) and (c) the presented topologies (strip and patch) allow for synthetizing low impedances only for unfeasible values of their structural parameters (i.e., g and w) in the THz range. For

strips narrower than 3 μm. As it will appear clearly in Section 2.6, this value could present some issues during the photolithography fabrication process [15]. Moreover, the high variation of the reactance with respect to the angle of incidence would produce very different radiating

However, several HMS geometries with limited spatial dispersion for both TE and TM polarizations have been proposed in the literature [35, 39], even for the development of highly directive FPC-LWAs [15]. Among them, the so-called fishnet geometry (see Figure 3(a)), i.e., a superposition between the patch and the strip layouts, will be investigated in Section 3. We note that the synthesis of HMSs with limited spatial dispersion is a key aspect for the development of reconfigurable LWAs. In fact, it allows for designing and optimizing the radiating

Figure 3. (a) Unit cells of three different type of HMSs: the strip grating, the array of patches, and the fishnet. (b) and (c) Surface reactance Xs vs. frequency in the range 0:5 < f < 1:5 [THz] for angles of incidence in the range 0<sup>∘</sup> ≤ θ ≤ 80<sup>∘</sup> with angular steps of 10<sup>∘</sup> : (b) strip grating for TM polarization and (c) array of patches for TE polarization (parameters in the text).

, an impedance of 20 Ω is reached with a strip grating with

), without affecting the

assuming ε<sup>r</sup> ¼ 2:3.

96 Metamaterials and Metasurfaces

example, at 1 THz and for <sup>θ</sup> <sup>¼</sup> <sup>0</sup><sup>∘</sup>

features as the beam is scanned from broadside to endfire.

impedance properties as the beam is steered.

properties of a reconfigurable FPC-LWA at broadside (i.e., <sup>θ</sup> <sup>¼</sup> <sup>0</sup><sup>∘</sup>

In a substrate-superstrate configuration, the PRS is represented by a cover layer made of a denser dielectric material with respect to the substrate, as depicted in Figure 1(b). These kinds of structures, commonly known as substrate-superstrate (SS) LWAs, have been extensively studied by N. G. Alexópoulos and D. R. Jackson [27, 40]. In those works [27, 40], they showed that the PRS (represented by the superstrate) would act as an almost perfectly reflecting layer, when the thickness h<sup>2</sup> of the superstrate (the cover layer) is chosen so that it is an odd multiple of one-quarter of the wavelength in the vertical direction. This corresponds to the condition

$$h\_2 = (2m - 1)\lambda/4 \left(\varepsilon\_{r2} - \sin^2 \theta\_0\right)^{-1/2}.\tag{9}$$

It is clear from Eqs. (4) and (9) that the design rules for optimizing radiation at broadside (θ<sup>0</sup> ¼ 0) with the fundamental (m ¼ 1) TE, TM leaky modes suggest a thickness of a half- and of a quarter-wavelength in their respective media for the substrate and the superstrate, respectively.

Another distinctive feature of SS-LWAs is that the directivity increases as the permittivity of the superstrate layer increases with respect to that of the substrate layer since the superstrate PRS acts as a more reflective surface.

#### 2.4.3. Multilayered dielectric structure

As shown in [27, 28], a very high step index (i.e., ffiffiffiffiffiffiffiffiffiffiffiffi εr2=ε<sup>r</sup> p ) is required for achieving relatively high directivities with a substrate-superstrate configurations. As high-permittivity materials usually exhibit also nonnegligible loss tangent, the radiation efficiency of a substrate-superstrate FPC-LWA is often limited by the relatively high dielectric losses of the dielectric layers. In this regard, a different way to increase directivity was proposed in [41]. There, the single superstrate is replaced by a periodic array of such superstrates, as shown in Figure 1(c). In this structure, the PRS consists of a stack of multiple alternating layers of high permittivity and low permittivity, whose thicknesses are chosen according to Eq. (9). As shown in [41], the directivity increases geometrically with the number of superstrate layers, and thus very directive beams may be obtained using modest values of superstrate permittivity, provided that several superstrate layers are used. From a circuital point of view, the gain enhancement phenomenon of this class of structures is readily explained by noticing that the 2n � 1 quarter-wavelength sections reduce the free-space impedance (seen as a load impedance of the transmission line segment looking upward from the substrate interface) of a factor ð Þ <sup>ε</sup>r2=ε<sup>r</sup> <sup>n</sup> .

As a last remark, it is worth noting that the multistack of alternating layers can be interpreted as a distributed Bragg reflector (DBR) meaning that it is a periodic electromagnetic bandgap (EBG) structure operating in a stopband. From this point of view, the substrate acts as a defect in the periodic EBG structure [42]. More interestingly, a leaky-wave explanation of the phenomenon is provided in [29].

#### 2.5. Reconfigurable FPC-LWAs

In most of microwave applications, the frequency scanning behavior typical of conventional FPC-LWAs is an undesirable effect, since the antenna radiating performance should not change within a small fractional bandwidth centered around the operating frequency. Conversely, it would be extremely beneficial, e.g., for SATCOM applications to dynamically change certain radiating properties (e.g., the pointing angle) at a fixed frequency. This would allow for avoiding solutions based on either heavy mechanical rotation systems or expensive 2-D phased arrays.

two types of excitations suggest to realize a quasi-resonant slot (thus modeling a half-wavelength dipole rather than a short dipole) for a twofold reason. With regard to free-space excitation, one has to consider that commercially available THz lenses focus on the energy over a spot size around 1 mm with a 2-D Gaussian profile [50]. As a result, a good energy coupling is obtained if one considers a slot of dimensions not less than 100 μm; such a dimension corresponds to λ=3 at the targeted frequency of 1 THz. With regard to guided wave excitation, one has to consider the cross section of commercial THz waveguides operating in the THz range. In the 900–1400 GHz band, the waveguide cross section is 200 μm � 100 μm [15], corresponding to 2λ=3 � λ=3 at

Terahertz Leaky-Wave Antennas Based on Metasurfaces and Tunable Materials

http://dx.doi.org/10.5772/intechopen.78939

99

In Section 2.4, we have seen that the PRS of an FPC-LWA may consist of a homogenized metasurface (HMS). As an HMS consists of a periodic lattice of subresonant elements, patterning of these elements on a metallic sheet is subject to the fabrication tolerances. Since the smallest details commonly allowed in a standard photolithographic process should be greater than 3 μm to avoid fabrication issues in low-cost, large-area production [15], the unit-cell period is set to p ¼ λ0=5 ¼ 60 μm. On the one hand, this value fulfills the homogenization limit [34] (p should preferably be smaller than λ0=4). On the other hand, the choice of p ¼ λ0=5 allows for designing details of the mask down to p=20 when designing the HMS at 1 THz.

The realization of a biasing scheme in a reconfigurable THz FPC-LWA based on graphene differs from that of an FPC-LWA based on nematic liquid crystals (NLCs). In an NLC-based FPC-LWA, a pair of conductive layers "sandwich" the NLC cell. However, a conductive layer introduces ohmic losses in the cavity; hence, it is important to use quasi-transparent polymers in the THz range. We therefore chose the PEDOT:PSS since it is a conductive film with a moderate conductivity in the THz range and an extremely thin profile [15]. Conversely, in a graphene-based FPC-LWA, a pair of electrodes is placed on top of the graphene sheet and on a conductive layer placed underneath of it. The conductive layer might be either a PEDOT:PSS [20], or a polysilicon layer [45]. In any case, it is important that the material between graphene and the conductive layer exhibits a high-voltage breakdown, so that graphene chemical potential can be raised at sufficiently high values to considerably change the conductive properties of the graphene sheet [20]. Examples of materials exhibiting a high-voltage breakdown are TiO2, Al2O3, and HfO2 [20]. As a final comment, we note that the thicknesses of the conductive layers and of the dielectric filling the gap with the graphene sheet are electrically small, and

can therefore safely be neglected in the transmission line model of an FPC-LWA.

A fishnet-like unit-cell can be obtained as a superposition of patches and strips of the same dimensions (see Figure 3(a)). Differently from array of patches and strip grating, this geometry

3. FPC-LWAs based on fishnet metasurfaces

1 THz.

2.6.3. THz fabrication tolerance

2.6.4. THz biasing scheme

3.1. Fishnet properties

The idea at the root of a reconfigurable FPC-LWA is instead simple and effective. In order to obtain the beam-steering capability in an FPC-LWA, it is necessary to change the phase constant of the relevant leaky mode excited in the cavity at a fixed frequency. This can be done, e.g., by dynamically changing the relevant geometrical or electromagnetic parameters of the materials, i.e., the thicknesses of the layers, or the relative permittivities and permeabilities (surface impedances for zero-thickness layers), respectively. While this possibility has extensively been addressed in the microwave range through, e.g., active impedances [43], or ferroelectric and ferromagnetic materials [44], few realizations exist in the THz range [20, 22, 45, 46].

In those works [20, 22, 45, 46], a tunable metasurface (e.g., graphene [20, 45, 46]) or a tunable material (e.g., nematic liquid crystals [22]) are employed as shown in Figures 1(d)–(f), respectively. The idea is to exploit the voltage dependence of either the surface conductivity of graphene (ambipolar field effect [47]) or the dielectric permittivity of a nematic liquid crystal (nematic phase transition [17]). Therefore, by changing the electromagnetic properties of the material, one can affect the leaky wavenumber (and, in turn, the radiating beam properties of an FPC-LWA) by applying a control voltage through a biasing scheme. The realization of such biasing schemes in the THz range is challenging and thus motivates the lack of reconfigurable THz FPC-LWA realizations. This and other aspects concerning the severe constraints posed by THz technology are discussed in the next Section 2.6.

#### 2.6. FPC-LWAs in the THz range

When a reconfigurable FPC-LWA has to be designed in the THz range, one has to consider several additional technological constraints: (i) the availability of low loss THz materials, (ii) the availability of efficient THz sources, (iii) the fabrication tolerances, and (iv) the complexity in the realization of a biasing scheme. These aspects are discussed below on a point-by-point basis.

#### 2.6.1. THz materials

An FPC-LWA consists of a grounded dielectric substrate of thickness 0:5λ= ffiffiffiffi εr p covered with a PRS. In order to limit the losses in the cavity, the dielectric material should preferably exhibit a low loss tangent in the THz range. Examples of such materials are Zeonor (ε<sup>r</sup> ¼ 2:3 and tan δ ¼ 0:001) and quartz (ε<sup>r</sup> ¼ 3:8 and tan δ ¼ 0:014). The latter, even if characterized by a loss tangent greater by an order of magnitude with respect to that of Zeonor, is particularly amenable for the realization of graphene-based structures. Indeed, graphene, when transferred on SiO2, shows limited phonon scattering [48], an effect that heavily contributes to decreasing the quality of graphene during its synthesis (see [49] for further details).

#### 2.6.2. THz sources

The FPC-LWAs investigated here are assumed to be excited with an HMD. This type of excitation can be realized by etching a subresonant slot in the ground plane. Such a slot can be either back-illuminated with a coherent THz source, or fed with a THz waveguide. These two types of excitations suggest to realize a quasi-resonant slot (thus modeling a half-wavelength dipole rather than a short dipole) for a twofold reason. With regard to free-space excitation, one has to consider that commercially available THz lenses focus on the energy over a spot size around 1 mm with a 2-D Gaussian profile [50]. As a result, a good energy coupling is obtained if one considers a slot of dimensions not less than 100 μm; such a dimension corresponds to λ=3 at the targeted frequency of 1 THz. With regard to guided wave excitation, one has to consider the cross section of commercial THz waveguides operating in the THz range. In the 900–1400 GHz band, the waveguide cross section is 200 μm � 100 μm [15], corresponding to 2λ=3 � λ=3 at 1 THz.

#### 2.6.3. THz fabrication tolerance

would be extremely beneficial, e.g., for SATCOM applications to dynamically change certain radiating properties (e.g., the pointing angle) at a fixed frequency. This would allow for avoiding solutions based on either heavy mechanical rotation systems or expensive 2-D phased arrays.

The idea at the root of a reconfigurable FPC-LWA is instead simple and effective. In order to obtain the beam-steering capability in an FPC-LWA, it is necessary to change the phase constant of the relevant leaky mode excited in the cavity at a fixed frequency. This can be done, e.g., by dynamically changing the relevant geometrical or electromagnetic parameters of the materials, i.e., the thicknesses of the layers, or the relative permittivities and permeabilities (surface impedances for zero-thickness layers), respectively. While this possibility has extensively been addressed in the microwave range through, e.g., active impedances [43], or ferroelectric and ferromagnetic materials [44], few realizations exist in the THz range [20, 22, 45, 46].

In those works [20, 22, 45, 46], a tunable metasurface (e.g., graphene [20, 45, 46]) or a tunable material (e.g., nematic liquid crystals [22]) are employed as shown in Figures 1(d)–(f), respectively. The idea is to exploit the voltage dependence of either the surface conductivity of graphene (ambipolar field effect [47]) or the dielectric permittivity of a nematic liquid crystal (nematic phase transition [17]). Therefore, by changing the electromagnetic properties of the material, one can affect the leaky wavenumber (and, in turn, the radiating beam properties of an FPC-LWA) by applying a control voltage through a biasing scheme. The realization of such biasing schemes in the THz range is challenging and thus motivates the lack of reconfigurable THz FPC-LWA realizations. This and other aspects concerning the severe constraints posed by

When a reconfigurable FPC-LWA has to be designed in the THz range, one has to consider several additional technological constraints: (i) the availability of low loss THz materials, (ii) the availability of efficient THz sources, (iii) the fabrication tolerances, and (iv) the complexity in the realization of a biasing scheme. These aspects are discussed below on a point-by-point basis.

PRS. In order to limit the losses in the cavity, the dielectric material should preferably exhibit a low loss tangent in the THz range. Examples of such materials are Zeonor (ε<sup>r</sup> ¼ 2:3 and tan δ ¼ 0:001) and quartz (ε<sup>r</sup> ¼ 3:8 and tan δ ¼ 0:014). The latter, even if characterized by a loss tangent greater by an order of magnitude with respect to that of Zeonor, is particularly amenable for the realization of graphene-based structures. Indeed, graphene, when transferred on SiO2, shows limited phonon scattering [48], an effect that heavily contributes to decreasing

The FPC-LWAs investigated here are assumed to be excited with an HMD. This type of excitation can be realized by etching a subresonant slot in the ground plane. Such a slot can be either back-illuminated with a coherent THz source, or fed with a THz waveguide. These

εr

p covered with a

An FPC-LWA consists of a grounded dielectric substrate of thickness 0:5λ= ffiffiffiffi

the quality of graphene during its synthesis (see [49] for further details).

THz technology are discussed in the next Section 2.6.

2.6. FPC-LWAs in the THz range

98 Metamaterials and Metasurfaces

2.6.1. THz materials

2.6.2. THz sources

In Section 2.4, we have seen that the PRS of an FPC-LWA may consist of a homogenized metasurface (HMS). As an HMS consists of a periodic lattice of subresonant elements, patterning of these elements on a metallic sheet is subject to the fabrication tolerances. Since the smallest details commonly allowed in a standard photolithographic process should be greater than 3 μm to avoid fabrication issues in low-cost, large-area production [15], the unit-cell period is set to p ¼ λ0=5 ¼ 60 μm. On the one hand, this value fulfills the homogenization limit [34] (p should preferably be smaller than λ0=4). On the other hand, the choice of p ¼ λ0=5 allows for designing details of the mask down to p=20 when designing the HMS at 1 THz.

#### 2.6.4. THz biasing scheme

The realization of a biasing scheme in a reconfigurable THz FPC-LWA based on graphene differs from that of an FPC-LWA based on nematic liquid crystals (NLCs). In an NLC-based FPC-LWA, a pair of conductive layers "sandwich" the NLC cell. However, a conductive layer introduces ohmic losses in the cavity; hence, it is important to use quasi-transparent polymers in the THz range. We therefore chose the PEDOT:PSS since it is a conductive film with a moderate conductivity in the THz range and an extremely thin profile [15]. Conversely, in a graphene-based FPC-LWA, a pair of electrodes is placed on top of the graphene sheet and on a conductive layer placed underneath of it. The conductive layer might be either a PEDOT:PSS [20], or a polysilicon layer [45]. In any case, it is important that the material between graphene and the conductive layer exhibits a high-voltage breakdown, so that graphene chemical potential can be raised at sufficiently high values to considerably change the conductive properties of the graphene sheet [20]. Examples of materials exhibiting a high-voltage breakdown are TiO2, Al2O3, and HfO2 [20]. As a final comment, we note that the thicknesses of the conductive layers and of the dielectric filling the gap with the graphene sheet are electrically small, and can therefore safely be neglected in the transmission line model of an FPC-LWA.

#### 3. FPC-LWAs based on fishnet metasurfaces

#### 3.1. Fishnet properties

A fishnet-like unit-cell can be obtained as a superposition of patches and strips of the same dimensions (see Figure 3(a)). Differently from array of patches and strip grating, this geometry exhibits almost negligible spatial dispersion [15]. Therefore, the surface impedance of a fishnetlike unit-cell is accurately represented by a purely imaginary scalar impedance value Zs ¼ jXs. As a matter of fact, this geometry has many other advantages with respect to the previous configurations [15]. As shown in [15], the fishnet design allows for spanning values of Xs from 100 Ω to few ohms in the same parameter space. Thus, highly directive FPC-LWAs can be designed using a suitable combination of g and w values. In fact, since the fishnet element has a two-valued space parameter, the designer has more degrees of freedom to synthesize the required values of the impedance. For example, a reactance as low as Xs ¼ 15 Ω can be synthesized by a fishnet unit-cell with w ¼ 0:4p, g ¼ 0:2p, and p ¼ λ0=5 ¼ 60 μm, corresponding to w ¼ 24 μm and g ¼ 12 μm, and thus fully compatible with a low-cost standard photolithographic fabrication process, as discussed in Section 2.6.

Therefore, the proposed fishnet element appears to be a valuable choice for the realization of the HMS since it allows for a flexible synthesis of the surface impedance. Moreover, its remarkably low spatial dispersion allows for a simple but accurate modeling, as will be shown in the next Section 3.2.

simulations: one of the ideal, lossless structure to determine αrad, and one of the lossy structures to determine α<sup>z</sup> and in turn αloss. With this method on hand, we obtained er ¼ 69% for L ¼ 20λ<sup>0</sup>

validated with full-wave simulations (blue circles) are reported over the (b) E-plane and (c) H-plane [15].

Figure 4 (a) Dispersion curves βz=k<sup>0</sup> (solid lines) and αz=k<sup>0</sup> vs. frequency (dashed lines) in the spectral range from 0:5 to 1:5 THz for TE (in red) and TM (in blue) polarizations. The inset highlights the behavior at the leaky cutoff, i.e., around f ¼ 1 THz. (b) and (c) Normalized radiation patterns calculated with the leaky-wave approach (black solid lines) and

Terahertz Leaky-Wave Antennas Based on Metasurfaces and Tunable Materials

http://dx.doi.org/10.5772/intechopen.78939

101

Graphene is a one-atom-thick layer of carbon atoms arranged in a honeycomb lattice with exceptional electronic transport properties [47]. From an antenna engineering point of view, one extremely interesting aspect of graphene theory is that due to its infinitesimal thickness, a graphene monolayer is adequately treated as a metasurface whose homogenized surface conductivity σ (neglecting nonlocal effects [51]) can be derived in scalar form by means of the Kubo formalism [47]. In the low THz range, i.e., for 0:3 ≤ f ≤ 3 THz and at room temperature, i.e., T ¼ 300 K, the complex-valued surface conductivity of graphene σ ¼ σ<sup>R</sup> � jσ<sup>J</sup> (where σ<sup>R</sup> and �σ<sup>J</sup> expressed in Siemens [S] are the conductance and the susceptance of graphene equivalent admittance, respec-

tively) is sufficiently well described by the following Drude-like expression [52, 53]:

<sup>e</sup> kBT

where ω ¼ 2πf is the angular frequency, �qe is the electron charge, kB is the Boltzmann constant, ℏ is the reduced Planck constant, τ is the relaxation time (which mainly determines graphene quality), and μ<sup>c</sup> is the chemical potential. In the following, a relaxation time τ ¼ 3 ps

Graphene quality strongly varies sample by sample, depending on several factors affecting the graphene synthesis. As a consequence, the relaxation time τ is never known a priori. The interested reader can refer to the recent detailed survey

<sup>τ</sup>ð Þ �<sup>1</sup> <sup>þ</sup> <sup>j</sup><sup>ω</sup> <sup>π</sup>ℏ<sup>2</sup> ln 2cosh <sup>μ</sup><sup>c</sup>

2kBT

, (10)

<sup>σ</sup> <sup>¼</sup> <sup>σ</sup><sup>R</sup> � <sup>j</sup>σ<sup>J</sup> <sup>¼</sup> <sup>2</sup>q<sup>2</sup>

(which leads to η<sup>r</sup> > 90%), and a gain greater than 26 dB.

4. FPC-LWAs based on graphene

4.1. Graphene properties

is assumed.<sup>4</sup>

proposed in [49] as well as discussions in [20, 21].

4

### 3.2. Fishnet-based design

The considered structure consists of a GDS made of Zeonor, covered with a fishnet-based HMS with parameters g ¼ 0:2p and w ¼ 0:4p, in order to synthesize an impedance of 15 Ω at 1 THz. While not spatially dispersive, the fishnet unit-cell is evidently frequency-dispersive, showing a quasi-linear frequency dependence [15], thus suggesting a dominant inductive behavior.

It is worth noting here that the negligible spatial dispersion of the fishnet allows for safely modeling the impedance Zsð Þf as a function of the frequency f only, without including the dependence on the wavenumber (see [15] for further details). This model of Zsð Þf has then been implemented in the equivalent circuit model depicted in Figure 2 (replacing ZPRS with Zsð Þf ), to perform the dispersion analysis of the structure.

Numerical results are reported in Figure 4(a) where the behaviors of β^ <sup>z</sup> (solid lines) and α^<sup>z</sup> (dashed lines) are shown in the frequency range 0:5–1:5 THz for both TE (in red) and TM (in blue) polarizations. As expected, the fishnet acts as a very high reflective PRS (due to the low impedance value) and the resulting leaky modes are slight perturbations of those of an equivalent dielectric-filled PPW. The very low value of α^<sup>z</sup> at the leaky cutoff β<sup>z</sup> ≃α<sup>z</sup> suggests a very narrow pencil beam. This result is confirmed by the radiation patterns obtained through leakywave theory (black solid lines) and CST full-wave simulations (blue circles), as shown in Figure 4(b) for the E-plane and in Figure 4(c) for the H-plane. The remarkable agreement between the results confirms that the TE-TM leaky-wave pair is dominant, and that the radiation efficiency η<sup>r</sup> is high.

In this regard, it is worth commenting on the overall expected efficiency er of this device to have a measure of the gain (G) rather than directivity (D0) (we recall that G ¼ erD0). As a matter of fact, the very low but finite loss tangent of Zeonor [22] contributes to raise the leakage rate. Hence, we have determined the amount of leakage due to radiation α^rad, and that due to absorption α^loss (where α^<sup>z</sup> ¼ α^loss þ α^radÞ by running two different numerical

Terahertz Leaky-Wave Antennas Based on Metasurfaces and Tunable Materials http://dx.doi.org/10.5772/intechopen.78939 101

Figure 4 (a) Dispersion curves βz=k<sup>0</sup> (solid lines) and αz=k<sup>0</sup> vs. frequency (dashed lines) in the spectral range from 0:5 to 1:5 THz for TE (in red) and TM (in blue) polarizations. The inset highlights the behavior at the leaky cutoff, i.e., around f ¼ 1 THz. (b) and (c) Normalized radiation patterns calculated with the leaky-wave approach (black solid lines) and validated with full-wave simulations (blue circles) are reported over the (b) E-plane and (c) H-plane [15].

simulations: one of the ideal, lossless structure to determine αrad, and one of the lossy structures to determine α<sup>z</sup> and in turn αloss. With this method on hand, we obtained er ¼ 69% for L ¼ 20λ<sup>0</sup> (which leads to η<sup>r</sup> > 90%), and a gain greater than 26 dB.

#### 4. FPC-LWAs based on graphene

#### 4.1. Graphene properties

exhibits almost negligible spatial dispersion [15]. Therefore, the surface impedance of a fishnetlike unit-cell is accurately represented by a purely imaginary scalar impedance value Zs ¼ jXs. As a matter of fact, this geometry has many other advantages with respect to the previous configurations [15]. As shown in [15], the fishnet design allows for spanning values of Xs from 100 Ω to few ohms in the same parameter space. Thus, highly directive FPC-LWAs can be designed using a suitable combination of g and w values. In fact, since the fishnet element has a two-valued space parameter, the designer has more degrees of freedom to synthesize the required values of the impedance. For example, a reactance as low as Xs ¼ 15 Ω can be synthesized by a fishnet unit-cell with w ¼ 0:4p, g ¼ 0:2p, and p ¼ λ0=5 ¼ 60 μm, corresponding to w ¼ 24 μm and g ¼ 12 μm, and thus fully compatible with a low-cost standard photolitho-

Therefore, the proposed fishnet element appears to be a valuable choice for the realization of the HMS since it allows for a flexible synthesis of the surface impedance. Moreover, its remarkably low spatial dispersion allows for a simple but accurate modeling, as will be shown in the next

The considered structure consists of a GDS made of Zeonor, covered with a fishnet-based HMS with parameters g ¼ 0:2p and w ¼ 0:4p, in order to synthesize an impedance of 15 Ω at 1 THz. While not spatially dispersive, the fishnet unit-cell is evidently frequency-dispersive, showing a quasi-linear frequency dependence [15], thus suggesting a dominant inductive behavior.

It is worth noting here that the negligible spatial dispersion of the fishnet allows for safely modeling the impedance Zsð Þf as a function of the frequency f only, without including the dependence on the wavenumber (see [15] for further details). This model of Zsð Þf has then been implemented in the equivalent circuit model depicted in Figure 2 (replacing ZPRS with Zsð Þf ),

(dashed lines) are shown in the frequency range 0:5–1:5 THz for both TE (in red) and TM (in blue) polarizations. As expected, the fishnet acts as a very high reflective PRS (due to the low impedance value) and the resulting leaky modes are slight perturbations of those of an equivalent dielectric-filled PPW. The very low value of α^<sup>z</sup> at the leaky cutoff β<sup>z</sup> ≃α<sup>z</sup> suggests a very narrow pencil beam. This result is confirmed by the radiation patterns obtained through leakywave theory (black solid lines) and CST full-wave simulations (blue circles), as shown in Figure 4(b) for the E-plane and in Figure 4(c) for the H-plane. The remarkable agreement between the results confirms that the TE-TM leaky-wave pair is dominant, and that the

In this regard, it is worth commenting on the overall expected efficiency er of this device to have a measure of the gain (G) rather than directivity (D0) (we recall that G ¼ erD0). As a matter of fact, the very low but finite loss tangent of Zeonor [22] contributes to raise the leakage rate. Hence, we have determined the amount of leakage due to radiation α^rad, and that due to absorption α^loss (where α^<sup>z</sup> ¼ α^loss þ α^radÞ by running two different numerical

<sup>z</sup> (solid lines) and α^<sup>z</sup>

graphic fabrication process, as discussed in Section 2.6.

to perform the dispersion analysis of the structure.

Numerical results are reported in Figure 4(a) where the behaviors of β^

Section 3.2.

3.2. Fishnet-based design

100 Metamaterials and Metasurfaces

radiation efficiency η<sup>r</sup> is high.

Graphene is a one-atom-thick layer of carbon atoms arranged in a honeycomb lattice with exceptional electronic transport properties [47]. From an antenna engineering point of view, one extremely interesting aspect of graphene theory is that due to its infinitesimal thickness, a graphene monolayer is adequately treated as a metasurface whose homogenized surface conductivity σ (neglecting nonlocal effects [51]) can be derived in scalar form by means of the Kubo formalism [47]. In the low THz range, i.e., for 0:3 ≤ f ≤ 3 THz and at room temperature, i.e., T ¼ 300 K, the complex-valued surface conductivity of graphene σ ¼ σ<sup>R</sup> � jσ<sup>J</sup> (where σ<sup>R</sup> and �σ<sup>J</sup> expressed in Siemens [S] are the conductance and the susceptance of graphene equivalent admittance, respectively) is sufficiently well described by the following Drude-like expression [52, 53]:

$$\delta\sigma = \sigma\_R - j\sigma\_\parallel = \frac{2q\_c^2 k\_B T}{(\pi^{-1} + j\omega)\pi\hbar^2} \ln\left[2\cosh\left(\frac{\mu\_c}{2k\_B T}\right)\right],\tag{10}$$

where ω ¼ 2πf is the angular frequency, �qe is the electron charge, kB is the Boltzmann constant, ℏ is the reduced Planck constant, τ is the relaxation time (which mainly determines graphene quality), and μ<sup>c</sup> is the chemical potential. In the following, a relaxation time τ ¼ 3 ps is assumed.<sup>4</sup>

<sup>4</sup> Graphene quality strongly varies sample by sample, depending on several factors affecting the graphene synthesis. As a consequence, the relaxation time τ is never known a priori. The interested reader can refer to the recent detailed survey proposed in [49] as well as discussions in [20, 21].

Clearly, σ is strongly affected by the values of μ<sup>c</sup> which in turn is related to the electrostatic bias E<sup>0</sup> through an integral equation [53]. As is known [20, 45, 53], the maximum absolute value of the chemical potential that can be obtained with an electrostatic field on the order of several V/nm (typically used in experiments5 ) is around 1 eV. It should be noted that such an electrostatic field is rather high and demands for the choice of suitable dielectric materials as discussed in Section 2.6.

Recent works [18, 20, 52, 53] have shown that in the THz range and for high-quality samples of graphene, the resistive part of graphene conductivity (σR) slightly increases as μ<sup>c</sup> increases, whereas its reactive part (σJ) considerably increases as μ<sup>c</sup> increases. In particular, for high values of μc, σ becomes mostly reactive so that graphene can be switched from a bad to a good conductor when μ<sup>c</sup> is raised in the range 0–1 eV. However, ohmic losses increase for high values of μ<sup>c</sup> [20]. Hence, biased graphene, even if of good quality (high values of τ), behaves as a good conductor with nonnegligible ohmic losses in the considered THz range.

Nevertheless, the possibility to considerably change the surface conductivity through the application of a bias voltage is one of the most exceptional aspects of graphene behavior since it allows for designing FPC-LWAs with beam-steering capabilities at fixed frequency. Such devices will be the object of the next Section 4.2.

The beam scanning capability through bias at fixed frequency is confirmed by the far-field radiation patterns (see colored dashed lines Figures 6(b) and (c)) obtained through CST fullwave simulations at f ¼ 0:922 THz. As shown, for values of the chemical potential spanning

Terahertz Leaky-Wave Antennas Based on Metasurfaces and Tunable Materials

http://dx.doi.org/10.5772/intechopen.78939

103

However, the beamwidths are not very narrow (and, in turn, the directivities are rather low) due to the relatively large attenuation constant. This means that a graphene sheet alone is not

To circumvent this issue, a slightly different configuration, also known as the graphene-based substrate-superstrate (GSS) antenna, has been proposed in the literature [19]. The idea underlying this novel configuration was the introduction of an additional cover layer on top of the GPW to increase the reflectivity of the overall PRS, now constituted by the combination of the graphene layer and the cover layer (see Figure 5(b)). Following the idea of a substratesuperstrate configuration, a quarter-wavelength dielectric material of high permittivity (namely, HfO2 with εr<sup>2</sup> ¼ 25, thus fixing the thickness to h<sup>2</sup> ≃14 μm) has been introduced on top of the GPW. However, the modal configuration of the tangential components of the electric field of an unperturbed substrate-superstrate structure (without the graphene sheet) has a minimum at the substrate-superstrate interface. As a consequence, a graphene sheet placed at the interface weakly interacts with the cavity and thus does not allow for achieving an effective reconfigurability of the pattern. This aspect poses the problem to find an optimal position for the graphene sheet within the cavity. As shown in [19], the position of the graphene sheet that maximizes the directivity in a GSS configuration is at x<sup>0</sup> ¼ 0:82h<sup>1</sup> (see Figure 5(b)) for an operating frequency of f ¼ 1:132 THz, i.e., in a middle way between a point of maximum interaction (i.e., the middle plane, x<sup>0</sup> ¼ 0:5h1) and a point of minimum interaction (i.e., the interface x<sup>0</sup> ¼ h1). The reason for this intermediate position is due to the nonnegligible ohmic losses of graphene, represented by the real part of its surface conductivity. Indeed, as the graphene sheet interacts more with the tangential electric field, the resonance effect increases (thus decreasing the leakage rate due to radiation α^rad), but also the contribution of the losses raises (thus increasing that part of the leakage rate accounting for the losses α^loss). As a result, the radiation efficiency and the directivity at broadside of a GSS strongly depend on the position of graphene inside the cavity. This has been quantitatively explained by the power analysis carried out in [20]. Unfortunately, the optimal position in terms of directivity at

, respectively.

from 1 to 0:5 eV, the GPW is able to steer the beam from broadside to 45<sup>∘</sup>

sufficiently reflective as a PRS, even when highly biased.

Figure 5. 2-D section and transmission line model of (a) a GPW and (b) a GSS [20].

#### 4.2. Graphene-based designs

In the simplest graphene-based FPC-LWA configuration, the PRS is replaced by a uniform (nonpatterned<sup>6</sup> ) graphene sheet (see Figure 5(a)), possibly biased as discussed in Section 2.6. We call this configuration a graphene planar waveguide (GPW) [18]. As in a conventional FPC-LWA antenna design [12], the thickness of the dielectric layer h is set according to Eq. (4). Considering a GPW filled with a SiO2 layer (ε<sup>r</sup> ¼ 3:8) and operated at f <sup>0</sup> ¼ 1 THz, it results h≃ 77 μm. The GPW is therefore similar to the fishnet-based FPC-LWA seen in Section 3, but the presence of the graphene sheet in place of the fishnet allows for achieving beam scanning at a fixed frequency by tuning the graphene surface conductivity through the biasing scheme.

As shown in [18], if one solves the dispersion equations (see Eq. (6)) by replacing YPRS with σ<sup>7</sup> , one finds that the phase constants of the fundamental TE, TM leaky mode pair considerably change as the chemical potential μ<sup>c</sup> is decreased from 1 to 0 eV. In particular, at f ¼ 0:922 THz (corresponding to the leaky cutoff condition of the TM mode of a GPW biased at μ<sup>c</sup> ¼ 1 eV [18]), the normalized phase constants of the leaky modes span values from about 0:2 to almost 1 for μ<sup>c</sup> varying from 1 to 0 eV (see blue and red dashed lines in Figure 6(a) for the dispersion curves of the TM and the TE leaky modes of the GPW, respectively). According to Eq. (3), this allows for steering the beam from broadside to endfire.

<sup>5</sup> In recent experiments, bias voltages Vb on the order of 100 V have been applied across a grounded dielectric slab covered with a graphene sheet [54]. In that case, a slab of quartz (ε<sup>r</sup> ¼ 3:8) of thickness t ¼ 300 μm was considered. Thus, for typical bias voltages in the range of 0 � 100 V, <sup>E</sup><sup>0</sup> <sup>¼</sup> Vg=<sup>t</sup> would be in the range 0 � 1 V/nm [54]. <sup>6</sup>

FPC-LWA designs based on patterned graphene sheets (namely, a 1-D array of capacitive strips of graphene) have been preliminarily investigated in [21].

<sup>7</sup> The same expression of σ can be used for both the TM and the TE polarization, since graphene exhibits negligible spatial dispersion, except for extremely slow waves such as surface plasmons [51], a case which is of no interest here.

Terahertz Leaky-Wave Antennas Based on Metasurfaces and Tunable Materials http://dx.doi.org/10.5772/intechopen.78939 103

Figure 5. 2-D section and transmission line model of (a) a GPW and (b) a GSS [20].

Clearly, σ is strongly affected by the values of μ<sup>c</sup> which in turn is related to the electrostatic bias E<sup>0</sup> through an integral equation [53]. As is known [20, 45, 53], the maximum absolute value of the chemical potential that can be obtained with an electrostatic field on the order of several V/nm

rather high and demands for the choice of suitable dielectric materials as discussed in Section 2.6. Recent works [18, 20, 52, 53] have shown that in the THz range and for high-quality samples of graphene, the resistive part of graphene conductivity (σR) slightly increases as μ<sup>c</sup> increases, whereas its reactive part (σJ) considerably increases as μ<sup>c</sup> increases. In particular, for high values of μc, σ becomes mostly reactive so that graphene can be switched from a bad to a good conductor when μ<sup>c</sup> is raised in the range 0–1 eV. However, ohmic losses increase for high values of μ<sup>c</sup> [20]. Hence, biased graphene, even if of good quality (high values of τ), behaves as

Nevertheless, the possibility to considerably change the surface conductivity through the application of a bias voltage is one of the most exceptional aspects of graphene behavior since it allows for designing FPC-LWAs with beam-steering capabilities at fixed frequency. Such

In the simplest graphene-based FPC-LWA configuration, the PRS is replaced by a uniform

We call this configuration a graphene planar waveguide (GPW) [18]. As in a conventional FPC-LWA antenna design [12], the thickness of the dielectric layer h is set according to Eq. (4). Considering a GPW filled with a SiO2 layer (ε<sup>r</sup> ¼ 3:8) and operated at f <sup>0</sup> ¼ 1 THz, it results h≃ 77 μm. The GPW is therefore similar to the fishnet-based FPC-LWA seen in Section 3, but the presence of the graphene sheet in place of the fishnet allows for achieving beam scanning at a fixed frequency by tuning the graphene surface conductivity through the biasing scheme. As shown in [18], if one solves the dispersion equations (see Eq. (6)) by replacing YPRS with σ<sup>7</sup>

one finds that the phase constants of the fundamental TE, TM leaky mode pair considerably change as the chemical potential μ<sup>c</sup> is decreased from 1 to 0 eV. In particular, at f ¼ 0:922 THz (corresponding to the leaky cutoff condition of the TM mode of a GPW biased at μ<sup>c</sup> ¼ 1 eV [18]), the normalized phase constants of the leaky modes span values from about 0:2 to almost 1 for μ<sup>c</sup> varying from 1 to 0 eV (see blue and red dashed lines in Figure 6(a) for the dispersion curves of the TM and the TE leaky modes of the GPW, respectively). According to Eq. (3), this allows for

In recent experiments, bias voltages Vb on the order of 100 V have been applied across a grounded dielectric slab covered with a graphene sheet [54]. In that case, a slab of quartz (ε<sup>r</sup> ¼ 3:8) of thickness t ¼ 300 μm was considered. Thus, for

FPC-LWA designs based on patterned graphene sheets (namely, a 1-D array of capacitive strips of graphene) have been

The same expression of σ can be used for both the TM and the TE polarization, since graphene exhibits negligible spatial

dispersion, except for extremely slow waves such as surface plasmons [51], a case which is of no interest here.

typical bias voltages in the range of 0 � 100 V, <sup>E</sup><sup>0</sup> <sup>¼</sup> Vg=<sup>t</sup> would be in the range 0 � 1 V/nm [54]. <sup>6</sup>

) graphene sheet (see Figure 5(a)), possibly biased as discussed in Section 2.6.

,

a good conductor with nonnegligible ohmic losses in the considered THz range.

) is around 1 eV. It should be noted that such an electrostatic field is

(typically used in experiments5

102 Metamaterials and Metasurfaces

4.2. Graphene-based designs

(nonpatterned<sup>6</sup>

5

7

devices will be the object of the next Section 4.2.

steering the beam from broadside to endfire.

preliminarily investigated in [21].

The beam scanning capability through bias at fixed frequency is confirmed by the far-field radiation patterns (see colored dashed lines Figures 6(b) and (c)) obtained through CST fullwave simulations at f ¼ 0:922 THz. As shown, for values of the chemical potential spanning from 1 to 0:5 eV, the GPW is able to steer the beam from broadside to 45<sup>∘</sup> , respectively. However, the beamwidths are not very narrow (and, in turn, the directivities are rather low) due to the relatively large attenuation constant. This means that a graphene sheet alone is not sufficiently reflective as a PRS, even when highly biased.

To circumvent this issue, a slightly different configuration, also known as the graphene-based substrate-superstrate (GSS) antenna, has been proposed in the literature [19]. The idea underlying this novel configuration was the introduction of an additional cover layer on top of the GPW to increase the reflectivity of the overall PRS, now constituted by the combination of the graphene layer and the cover layer (see Figure 5(b)). Following the idea of a substratesuperstrate configuration, a quarter-wavelength dielectric material of high permittivity (namely, HfO2 with εr<sup>2</sup> ¼ 25, thus fixing the thickness to h<sup>2</sup> ≃14 μm) has been introduced on top of the GPW. However, the modal configuration of the tangential components of the electric field of an unperturbed substrate-superstrate structure (without the graphene sheet) has a minimum at the substrate-superstrate interface. As a consequence, a graphene sheet placed at the interface weakly interacts with the cavity and thus does not allow for achieving an effective reconfigurability of the pattern. This aspect poses the problem to find an optimal position for the graphene sheet within the cavity. As shown in [19], the position of the graphene sheet that maximizes the directivity in a GSS configuration is at x<sup>0</sup> ¼ 0:82h<sup>1</sup> (see Figure 5(b)) for an operating frequency of f ¼ 1:132 THz, i.e., in a middle way between a point of maximum interaction (i.e., the middle plane, x<sup>0</sup> ¼ 0:5h1) and a point of minimum interaction (i.e., the interface x<sup>0</sup> ¼ h1). The reason for this intermediate position is due to the nonnegligible ohmic losses of graphene, represented by the real part of its surface conductivity. Indeed, as the graphene sheet interacts more with the tangential electric field, the resonance effect increases (thus decreasing the leakage rate due to radiation α^rad), but also the contribution of the losses raises (thus increasing that part of the leakage rate accounting for the losses α^loss). As a result, the radiation efficiency and the directivity at broadside of a GSS strongly depend on the position of graphene inside the cavity. This has been quantitatively explained by the power analysis carried out in [20]. Unfortunately, the optimal position in terms of directivity at

NLCs sandwiched between two substrates covered by electrodes. It is known that [17] if an alignment film is deposited on the substrates and no bias is applied (Vb ¼ 0 V, unbiased state), the optical axis of the LC molecules is parallel to the substrate plane, in the direction given by the alignment film (horizontal axis), whereas, when a sufficiently large driving voltage (Vb ¼ V∞, biased state) is applied across the LC, the optical axis is fully tilted and parallel to the applied electric field direction (vertical axis) [55]. While the maximum achievable tuning range can be inferred from the knowledge of these two limiting states, the voltage-dependent tunable properties of LCs require a rigorous study of the LC dynamics. This task can be performed by a numerical implementation of the Q-tensor theory of liquid crystals, a convenient theoretical formulation for the study of the LC orientation in confined geometries [56]. Among the different kinds of NLCs, we consider here the nematic mixture 1825 because of its high birefringence at THz frequencies [56]. The dielectric properties of such a material are

an NLC in its unbiased and fully biased states reduces to a diagonal matrix, whose diagonal elements are the ordinary ε<sup>o</sup> and extraordinary ε<sup>e</sup> relative permittivities of the NLC. For intermediate states, the NLC permittivity tensor is no longer a diagonal matrix. However, a thorough analysis of the Q-tensor has revealed that for an NLC aligned along the horizontal z-axis (see the reference frame in Figure 7), when a low-driving voltage is applied [56] across a fishnet unit-cell, the optical axis of the NLC tilts in the xz-plane exhibiting a negligible rotation over both the yz-

where εyy ¼ εxxð Þ¼ 0 εzzð Þ¼ V<sup>∞</sup> ε<sup>o</sup> ≃ 2:42 and εxxð Þ¼ V<sup>∞</sup> εzzð Þ¼ 0 ε<sup>e</sup> ≃ 3:76 for a z-oriented NLC in the THz range [56], respectively. It is worth noting here that the assumption of uniaxial crystal for the NLC layer allows for easily describing its behavior with a relatively simple equivalent circuit model [57]. On the contrary, when the off-diagonal components are nonnegligible, more complicated equivalent networks are needed [57, 58], also accounting for the dielectric tensor

The proposed device (see Figure 7(a)) consists of a multistack of alternating layers of thin alumina layers and of NLCs placed above a GDS. The choice of Zeonor for the substrate layer has been motivated by the index matching between its relative permittivity εr<sup>1</sup> ¼ 2:3 and the ordinary relative permittivity ε<sup>o</sup> ¼ 2:42 of the NLC layer in the THz range [56], as it is required

As we have seen in Section 2.4, the alternation of high- and low-permittivity layers, with thicknesses fixed at odd multiples of a quarter wavelength in their respective media, allows for obtaining a narrow beam radiating at broadside [29]. In the proposed device, the innovating feature is represented by the possibility of exploiting the tunable properties of the NLCs [56], here representing the low-permittivity layer. In particular, the application of a common control signal to the NLC layers allows for changing their dielectric properties, thus achieving beam-steering capability at a fixed frequency, similar to the one obtained with the graphene-

A more accurate model should take into account the non-zero values of the εxz and εzx components.

. As a uniaxial crystal, the permittivity tensor of

Terahertz Leaky-Wave Antennas Based on Metasurfaces and Tunable Materials

http://dx.doi.org/10.5772/intechopen.78939

105

, the NLC can still be locally modeled as a uniaxial

ð Þ Vb is a diagonal matrix diag εxxð Þ Vb ; εyy; εzzð Þ Vb

described by a complex permittivity tensor ε<sup>r</sup>

and the xy-plane. Thus, at first approximation8

to properly enhance the resonance condition in an FPC [29].

crystal whose complex permittivity tensor ε<sup>r</sup>

spatial distribution [55].

5.2. NLC-based designs

based FPC-LWAs shown in Section 4.

8

Figure 6 (a) Normalized phase constants and attenuation constants of the fundamental TM (in blue) and TE (in red) leaky modes of a GPW (dashed lines) at f ¼ 0:922 THz and of a GSS (solid lines) with graphene placed at the optimum position x<sup>0</sup> ¼ 0:82h<sup>1</sup> at a fixed frequency f ¼ 1:132 THz, as a function of the chemical potential in the range 1 > μc > 0 eV. (b) H-plane and (c) E-plane radiation patterns, normalized to the overall maximum (achieved at broadside), vs. elevation angle θ for the GSS antenna (solid lines) and for the GPW (dashed lines). The scanning behavior at a fixed frequency (f ¼ 1:132 THz for the GSS and f ¼ 0:92 THz for the GPW) is shown for four theoretical pointing angles θ<sup>0</sup> ¼ 0, 15, 30, 45<sup>∘</sup> . The chemical potentials for the GPW and the GSS are reported in the legend [19].

broadside (x<sup>0</sup> ¼ 0:82h1) does not lead to the best configuration in terms of radiation efficiency (x<sup>0</sup> ¼ h<sup>1</sup> or x<sup>0</sup> ¼ 0). However, the radiation performance improvement in terms of directivity of a GSS with respect to a GPW is paid at the expense of a moderate reduction of radiation efficiency. To give some numbers, a GPW shows D<sup>0</sup> ≃14 dB and er ≃ 69%, whereas an optimized GSS shows D<sup>0</sup> ≃18:5 dB and er ≃ 56%.

The dispersion analysis of a GSS and its radiating performance are reported in Figure 6(a) (see blue and red solid lines for the dispersion curves of the TM and the TE leaky modes of the GPW, respectively) and Figures 6(b) and (c) (see colored solid lines). As manifested by comparing the dashed and solid colored lines in Figures 6(b) and (c), the GSS shows considerably narrow beamwidths (thus higher directivities) over the entire beam scanning range. However, it should be noted that the GSS achieves the same angular range of a GPW over a wider range of variation of the chemical potential. This aspect reveals the slightly weaker degree of reconfigurability of the GSS with respect to the GPW.
