# Power balance and efficiency of metasurface antennas

M. Bodehou, D. Gonzalez-Ovejero, C. Craeye, S. Maci, I. Huynen, E. Martini

► To cite this version:

M. Bodehou, D. Gonzalez-Ovejero, C. Craeye, S. Maci, I. Huynen, et al., Power balance and efficiency of metasurface antennas. Scientific Reports, 2020, 10 (1), pp.17508. 10.1038/s41598-020-74674-w. hal-03002110

HAL Id: hal-03002110

https://hal.science/hal-03002110

Submitted on 6 Jan 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# Power balance and efficiency of metasurface antennas

Modeste Bodehou<fnref n="1" />✉, David González-Ovejero², Christophe Craeye<fnref n="1" />, Stefano Maci³, Isabelle Huynen<fnref n="1" /> & Enrica Martini³

This paper presents two methods for the efficient evaluation of the power balance in circular metasurface (MTS) antennas implementing arbitrary modulated surface impedances on a grounded dielectric slab. Both methods assume the surface current in the homogenized MTS to be known. The first technique relies on the surface current expansion with Fourier-Bessel basis functions (FBBF) and proceeds by integration of the Poynting vector on a closed surface. The second method is based on the evaluation of the residue of the electric field spectrum at the surface-wave (SW) pole, and is demonstrated by using a current expansion in Gaussian ring basis functions (GRBF). The surface current expansions can be directly obtained either by analyzing the antenna with a Method of Moments (MoM) tool for homogenized MTSs based on FBBF or GRBF, or derived by a projection process. From there, the power contributions, namely the total power delivered by the feed, the radiated power, the SW power, and the Ohmic power losses in the dielectric are computed. Several efficiency metrics are presented and discussed: tapering efficiency, conversion efficiency, loss factor, and diffraction factor. Since the MTS apertures at hand are leaky-wave (LW) antennas, the designer must find a compromise between the aperture efficiency and the conversion efficiency. This requires accurate and fast computational techniques for the efficiency. The present paper demonstrates for the first time that the efficiency of MTS antenna devices can be accurately evaluated in a few minutes. The compromise that should be made during the design process between the tapering efficiency and the conversion efficiency is highlighted. The impact on the efficiency of isotropic versus anisotropic MTS, uniform versus non-uniform modulation index, is analyzed. An excellent agreement is obtained between both approaches, commercial software, and experimental data.

Among different types of metasurfaces (MTSs)<fnref n="1" />⁻⁶, modulated MTS antennas are an emerging class of radiating apertures that exploit the modulation of a surface impedance in such a way that an excited cylindrical surface-wave (SW) is transformed into leaky-waves (LW)⁷⁻¹⁶. The SW is in general launched by a transverse magnetic (TM) feed, such as a simple monopole placed at the center of the MTS. The theoretical premises of this technology go back to the work of Oliner and Hessel¹⁷, in which a rigorous treatment of SW propagation over sinusoidally modulated impenetrable reactive boundary conditions was provided. Over the past decade, this work has inspired many researchers in the antenna community and has led to practical devices with high performance¹¹⁻¹⁴,¹⁷,¹⁸. For example, the introduction of anisotropy in the surface reactance in combination with a variation of the modulation depth and the local period have been used to synthesize a wide range of aperture fields with excellent beam polarization purity and tapering efficiency¹⁸⁻²⁰. In addition, very general reactance modulation (i.e. not necessarily locally sinusoidal) based on the solution of the electric field integral equation (EFIE) has been proposed recently to improve the beam shaping capability²¹ as well as multibeam operation²²,²³. After the synthesis of the required impedance modulation, the surface impedance is usually implemented using sub-wavelength patches, printed on a grounded dielectric slab²⁴. A progressive change of the parameters describing the patches (orientation, size, etc) allows one to mimic the required impedance on a Cartesian regular lattice. This implementation can be carried out quite accurately as can be seen from the comparison of the surface reactance and the physical MTS structure simulations²⁵,²⁶. This means that the antenna performance can be accurately predicted at the homogenized (surface impedance) level, which is much more computationally efficient, as demonstrated by the EFIE techniques in²⁶,²⁷. Those techniques provide an excellent computational performance through the discretization of the EFIE involving the grounded substrate Green's function²⁸ kernel, and using appropriate entire-domain basis functions. However, the usage of the substrate Green's function assumes the substrate to

<footnote n="1">¹ICTEAM Institute, Université catholique de Louvain, Place du Levant 3, 1348 Louvain-la-Neuve, Belgium. ²Univ Rennes, CNRS, Institut d'Electronique et de Télécommunications de Rennes (IETR), UMR 6164, 35000 Rennes, France. ³Department of Information Engineering and Mathematics, University of Siena, 53100 Siena, Italy. ✉email: modeste.bodehou@uclouvain.be</footnote>
<figure><img src="image_3.png" /><figcaption>Figure 1. Geometry for the metasurface antenna.</figcaption></figure>

be infinite, and therefore does not properly deal with possible reflection/diffraction effects that stem from the substrate truncation, when the SW is not sufficiently attenuated before reaching the substrate boundary. Dif- fraction and reflection occurring at the substrate's edges may produce sidelobes which are not predicted by the simulation models based on the infinite substrate Green's function when a significant power level reaches the MTS rim. Ohmic power losses may also be non negligible at high frequencies and should be incorporated in the analysis. It is therefore relevant to have a clear understanding, during the design process, of the power balance in MTS antennas.

In²⁹, the efficiency of MTS antennas has been treated with the flat optics (FO) formalism³⁰. However, this approach, although providing closed-form formulas, is limited to locally sinusoidal reactance modulation. The present paper aims at providing two efficient semi-analytical methods for the rapid and accurate evaluation of the efficiency of MTS antennas. The proposed methods are not necessarily limited to locally sinusoidal modulations. That means, the methods can handle arbitrary surface impedance profiles i.e. based on any synthesis method⁷⁻²³. The derivation will rely on the mathematical formalisms introduced in²⁶,²⁷, where the MTS is modeled as a sheet (penetrable) impedance on top of the substrate. This modeling allows one to finely take into account the substrate dispersion, and therefore provides more accurate results in comparison with the opaque (impenetrable) imped- ance modeling²⁵. In²⁷, the current distribution has been expanded into Fourier-bessel basis functions (FBBFs), while²⁶ uses Gaussian ring basis functions (GRBFs). The choice of those bases in the present paper, besides the fact they provide a fast MoM solution, is justified by their complementary properties. Along the radial coordinate, while FBBFs exhibit selectivity in the spectral domain, GRBFs are selective in space domain. Therefore, they are complementary in describing the fields features. In this paper, FBBFs are used to compute the total power delivered by the feed on the basis of the Poynting theorem. GRBFs current are exploited to compute the SW power flowing beyond the MTS using the residue theorem. Therefore, developing both approaches provides a good cross-validation methodology.

## Results

**Power contributions.** Illustration and methodology. Let us consider a circular-domain MTS printed on a grounded dielectric slab, as depicted in Fig. 1. The MTS is assumed in this paper to be fed at its center ($ρ = 0$) with a vertical elementary dipole. Nevertheless, the methods presented below are also applicable for an arbitrary excitation, assuming the excitation fields are a priori known.

Two different power descriptions are adopted depending on whether the substrate is truncated (Fig. 2a) or infinite (Fig. 2b), the latter occurring in the full wave analysis with the Green's function of the infinite slab. Although the schematization is approximate, it helps to understand the power balance. When the slab is trun- cated, the total power $P_{tot}$ delivered by the feed can be calculated by integrating the flux of the real part of the Poynting vector through the dashed red surface in Fig. 2a. This can be subdivided in the integration over the aperture $S_1$ and the one over the lateral rim $S_2$. We denote these two contributions as $P_{MTS}$ and $P_{rim}$, respectively. The summation of these two contributions gives the radiated power $P_{vis}^{trunc}$, namely the power associated with the visible contribution of the spectrum (far-field). We note that, although $P_{rim}$ is the main contribution associated with the diffracted field, it is not the only one, being also present the asymptotic end point contribution from $S_1$. If losses are present, the total power delivered by the feed can be written as $P_{tot} = P_{vis}^{trunc} + P_{Ω}$, where $P_{Ω}$ is the power lost in the dielectric and in the metal.

When the grounded substrate is infinite, it is convenient to use a surface of integration like the one in Fig. 2b, which cuts the dielectric substrate. Now, one more term occurs in the power balance, namely the power that is transported by the SW excited at the rim of the printed MTS. We call this term $P_{diff}^{sw}$; it is not visible in the far- field and therefore, the power balance is written as

$$P_{tot} = P_{vis} + P_{Ω} + P_{diff}^{sw}. \qquad (1)$$

Since $P_{tot}$ and $P_{Ω}$ are almost the same in the two cases, we may say that $P_{vis}^{trunc} \approx P_{vis} + P_{diff}^{sw}$. The contribution $P_{diff}^{sw}$ is in general undesired and leads to a limited conversion efficiency of the antenna. It is worth noting that
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Illustration of the power contributions. (a) Finite MTS over a finite grounded substrate. (b) Finite MTS over an infinite grounded substrate. The MTS is fed with a coax terminated by a circular patch.</figcaption></figure>

we do not distinguish here the visible (radiated) field from the surface, namely converted from the SW in LW, and the one radiated directly from the feed. In other terms, we are not interested here in estimating the “feed efficiency”, which is treated in$^{29}$.

The first goal of the present paper consists in evaluating the power contributions, i.e., estimating $P_{vis}$, $P_{diff}^{sw}$, $P_{tot}$, and $P_{\Omega}$ for a given MTS antenna and excitation. The first method, based on FBBFs, computes $P_{vis}$, $P_{\Omega}$, $P_{tot}$, then deduces $P_{diff}^{sw} = P_{tot} - P_{vis} - P_{\Omega}$. The second method, based on GRBFs, computes $P_{vis}$, $P_{diff}^{sw}$, and $P_{\Omega}$, and then derives $P_{tot} = P_{vis} + P_{\Omega} + P_{diff}^{sw}$. This means that the two methods are not only different from the used set of basis functions, but also in term of the physical power contributions that are directly estimated.

In the following, it is assumed that each component of the surface current **J** in a Cartesian unit vector system ($\hat{x}, \hat{y}$) is known in a Fourier-Bessel basis$^{27}$ or Gaussian ring basis$^{26}$

$$J(\rho, \phi) = \sum_{n=-N}^{N} \sum_{m=1}^{M} i_{mn}^{x} R_{m,n}^{x}(\rho, \phi) \hat{x} + i_{mn}^{y} R_{m,n}^{y}(\rho, \phi) \hat{y}, \quad (2)$$

where $R_{m,n}$ corresponds to a given FBBF or GRBF, and $\rho$ and $\phi$ are the radial and azimuthal coordinates respectively. Those classes of basis functions are defined over the whole circular domain and have been used in$^{26,27}$ for the analysis of the directivity of MTS antennas printed on an infinite grounded slab. One should note that both families of functions are defined in closed-form in spatial domain and admit closed-form Fourier transforms.

**Radiated (visible) power computation.** The power radiated by the MTS is derived by integrating the Poynting vector over an infinitely extended disk just above the MTS layer (x-y plane). The spectral domain electric field is denoted as $\tilde{E}(k_x, k_y)$, where $(k_x, k_y)$ are the spectral Cartesian coordinates. Using the Parseval theorem, this power can be evaluated in spectral domain as:

$$P_{vis} = \frac{1}{8\pi^2} \int_{0}^{2\pi} \int_{0}^{k_0} ||\tilde{\mathbf{E}}||^2 \frac{k_z}{k_0 \eta_0} k_{\rho} dk_{\rho} d\alpha, \quad (3)$$

where $(k_{\rho}, \alpha)$ are the spectral cylindrical coordinates, $k_z = \sqrt{k_0^2 - k_{\rho}^2}$, with $k_0$ being the free-space wavenumber, and $\eta_0$ is the free-space impedance. $||\tilde{\mathbf{E}}||$ is the norm of the spectral aperture electric field $\tilde{\mathbf{E}}$. The latter is the sum of the electric field radiated by the current induced on the MTS and that of the feed, namely:

$$\tilde{\mathbf{E}} = \tilde{\mathbf{G}}^{EJ} \cdot \tilde{\mathbf{J}} + \tilde{\mathbf{E}}_{i}, \quad (4)$$

where $\tilde{\mathbf{G}}^{EJ}$ is the appropriate dyadic spectral Green's function of the grounded slab and $\tilde{\mathbf{E}}_{i}$ is the spectrum of the forced total electric field over the surface in absence of MTS. The integral in (3) can be efficiently calculated using the closed-form spectrum of the GRBFs [26, Expression (42)] or FBBFs [27, Expression (4)].

**Surface-wave power computation.** This contribution is relevant to the power of the SW that proceeds beyond the MTS. In the infinite-slab model this power is trapped in the slab, and does not contribute to radiation. However, in practical realizations it is reflected or diffracted at the slab truncation, and therefore contributes to the radiation pattern in an uncontrolled manner. For this reason, it is important to minimize this trapped SW power, during the antenna design. In the following, a closed-form expression is derived to evaluate this contribution as a function of previously computed currents. To this end, suppose that the z-component of the spectrum of the aperture electric field is a regular function that can be written as a Fourier series versus spectral angle $\alpha$:

$$\tilde{E}_z(k_\rho, \alpha) = \sum_n \tilde{E}_{zn}(k_\rho) e^{-j n \alpha}. \quad (5)$$

This form automatically comes out from (4) if GRBFs or FBBFs are used for the MTS current representation. Furthermore, the feed only contributes to the 0-indexed term, due to its axial symmetry. The spatial field is
obtained by taking the inverse Fourier transform of the expression in (5). As summarized in the Appendix, for $\rho > a$, where a is the radius of the MTS, this can be divided into two contributions: the radiated field and the SW field. The latter is given by the contribution of the residue at the SW pole of the bare grounded slab, $\beta_0^{sw}$, and reads (see proof in the supplementary material)

$$E_z^{SW} = -\frac{j}{2}e^{-j\frac{\pi}{4}}\sqrt{\frac{2\beta_0^{sw}}{\pi\rho}}e^{-j\beta_0^{sw}\rho}\sum_n R_n e^{-jn\phi}, \quad (6)$$

where we have assumed $\beta_0^{sw} \rho \gg 1$ and

$$R_n = \lim_{k_\rho \to \beta_0^{sw}} (k_\rho - \beta_0^{sw}) \tilde{E}_{zn} = \tilde{J}_n^{TM}(\beta_0^{sw}) R_{GF} + \delta_{n,0} R_{feed}. \quad (7)$$

In (7), $\tilde{J}_n^{TM}(\beta_0^{sw}) = \int_0^{2\pi} \tilde{J}^{TM}(\beta_0^{sw}, \alpha) e^{j\alpha\phi} d\alpha$, $R_{GF}$ is the residue of the Green's function relating the TM currents on the MTS to the vertical electric field on top of the slab, and $R_{feed}$ is the residue of the field radiated by the feed on the bare slab. From (6), and taking into account the local structure of the SW electromagnetic field, the real part of the relevant Poynting vector is derived, and integrated over a cylindrical surface surrounding the radiating aperture (extending vertically from the ground plane to infinity), to obtain the following expression for the residual SW power

$$P_{diff}^{sw} = \frac{\omega\epsilon_0}{4\epsilon_r} \left[ \frac{\tan(k_{zd}h)}{k_{zd}} + \frac{h}{\cos^2(k_{zd}h)} + \frac{\epsilon_r}{\alpha_{za}} \right] \sum_n |R_n|^2, \quad (8)$$

where $\epsilon_0$ is the free-space permittivity, $\epsilon_r$ is the relative permittivity of the substrate, $\omega = 2\pi f$ is the angular frequency, $h$ is the slab thickness, $k_{zd}^2 = k_0^2\epsilon_r - (\beta_0^{sw})^2$ and $\alpha_{za}^2 = (\beta_0^{sw})^2 - k_0^2$. Notice that the evaluation of (8) only requires the knowledge of the residues in (7), which can be calculated for any given current spectrum. The general expression can be also specialized to the particular case of a given current representation. For instance, if GRBFs are used to represent MTS currents, it results

$$\tilde{J}_n^{TM}(\beta_0^{sw}) = \pi j^n \sum_{m=1}^{M} \{ [j i_{m,n+1}^x + i_{m,n+1}^y] \Phi_{m,n+1}(\beta_0^{sw}) + [i_{m,n-1}^y - i_{m,n-1}^x] \Phi_{m,n-1}(\beta_0^{sw}) \}, \quad (9)$$

where the $\Phi$ functions are reported in [26, Expression (3)].

**Power delivered by the feed.** Referring to the Poynting theorem, the total power delivered by a matched feeding source in absence of losses can be computed after integrating the Poynting vector along a closed surface comprising the source. We propose to use the closed surface represented in dotted red lines in Fig. 2b. This surface includes the MTS aperture ($S_1$), the rim of the antenna ($S_2$), and the ground plane. Since the ground plane is assumed to be a perfect electric conductor, its contribution is zero. Therefore, the power delivered by the feeder can be computed as:

$$P_{tot} = \frac{1}{2} \mathrm{Re} \left\{ \iint_{S_1} \mathbf{E}_t \times \mathbf{H}_t^* \cdot \hat{\mathbf{n}} \, dS + \iint_{S_2} \mathbf{E} \times \mathbf{H}^* \cdot \hat{\mathbf{\rho}} \, dS + \iiint_V \sigma ||\mathbf{E}||^2 dV \right\} = P_{MTS} + P_{rim} + P_{\Omega}, \quad (10)$$

where $\mathbf{H}$ is the magnetic field, $t, \ast$, and $\mathrm{Re}$ stand respectively for the tangential part, the complex conjugate operator, and the real part. $V$ is the volume enclosed by the surface $S_1 \cup S_2$ and the ground plane. Finally, $\sigma$ is the conductivity of the dielectric substrate. The first term ($P_{MTS}$) is referred to as the MTS contribution, the second one ($P_{rim}$) as the rim contribution, and the last term ($P_{\Omega}$) is referred to as the Ohmic power losses.

**MTS contribution.** The MTS contribution is explicitly rewritten as:

$$P_{MTS} = \frac{1}{2} \mathrm{Re} \left\{ \iint_{S_1} E_x \cdot H_y^* - E_y \cdot H_x^* dS \right\}, \quad (11)$$

where $E_{x,y}$ and $H_{x,y}$ are the Cartesian components of the electric and magnetic field respectively. The aperture tangential electric field tested with the complex conjugate of the FBBFs can be computed as:

$$[\mathbf{E}_t] = [\mathbf{Z}_{IBC}] [\mathbf{I}], \quad (12)$$

where $[\mathbf{Z}_{IBC}]$ and $[\mathbf{I}]$ are respectively the impedance boundary condition (IBC) matrix and the surface current vector. The IBC matrix is calculated with the MoM procedure described in [27]. Since FBBFs are orthogonal, the bases coefficients $e_{mn}$ of $[\mathbf{E}_t]$ into FBBFs can be obtained from the tested fields $e_{mn}^{\mathrm{tested}}$ as: $e_{mn} = e_{mn}^{\mathrm{tested}}/K(m, n)$, where $K(m, n)$ is a normalization factor [27, Equ. 35]. The aperture tangential magnetic field tested with FBBFs is computed as

$$[\mathbf{H}_t] = [\mathbf{Z}_G^H] [\mathbf{I}] - [\mathbf{V}^H], \quad (13)$$

where $[\mathbf{Z}_G^H]$ is the substrate matrix relating the surface current to the aperture magnetic field. $[\mathbf{V}^H]$ is the excitation magnetic field at the MTS plane, tested with the complex conjugate of the FBBFs. $[\mathbf{Z}_G^H]$ and $[\mathbf{V}^H]$ are computed using a procedure similar to that used for the computation of $[\mathbf{Z}_G]$ and $[\mathbf{V}]$ in [27].
Similarly to the electric field, the magnetic field coefficients into FBBFs are obtained from the tested fields through renormalization. The aperture field is therefore written in the FBBF basis as

$$E_t(\rho, \phi) = \sum_{n=-N}^{N} \sum_{m=1}^{M} e_{mn}^x R_{m,n}^x(\rho, \phi) \hat{x} + e_{mn}^y R_{m,n}^y(\rho, \phi) \hat{y} \quad (14)$$

$$H_t(\rho, \phi) = \sum_{n=-N}^{N} \sum_{m=1}^{M} h_{mn}^x R_{m,n}^x(\rho, \phi) \hat{x} + h_{mn}^y R_{m,n}^y(\rho, \phi) \hat{y}. \quad (15)$$

Inserting expressions (14) and (15) into (11), and using the orthogonality relation of the FBBFs (see [27, Appendix A]), the MTS contribution to the delivered power is given by

$$P_{MTS} = \frac{1}{2} \sum_{n=-N}^{N} \sum_{m=1}^{M} K(m, n) \text{Re} \{ e_{mn}^x h_{mn}^{y*} - e_{mn}^y h_{mn}^{x*} \}. \quad (16)$$

This means that the calculation of the MTS contribution to the power delivered by the feed can be carried out extremely fast (in practice in less than 1 min) once the currents on the MTS are known in term of FBBFs (e.g. using the MoM).

**Rim contribution.** The rim contribution in (10) can be developed as

$$P_{rim} = \frac{a}{2} \text{Re} \left\{ \int_{z=-h}^{0} \int_{\phi=0}^{2\pi} E_{\phi}(\rho=a, \phi, z) H_{z}^{*}(\rho=a, \phi, z) - E_{z}(\rho=a, \phi, z) H_{\phi}^{*}(\rho=a, \phi, z) d\phi \, dz \right\}. \quad (17)$$

Since the mode is quasi-TM over a ground plane, it is reasonable to neglect the first term. In the following, we will hence detail the calculation of the second term only; bearing in mind that the first term can be computed similarly.

$$P_{rim} \approx -\frac{a}{2} \text{Re} \left\{ \int_{z=-h}^{0} \int_{\phi=0}^{2\pi} E_z(\rho=a) H_{\phi}^{*}(\rho=a) d\phi \, dz \right\}, \quad (18)$$

where the arguments $\phi$ and $z$ have been omitted. Expression (18) can be rewritten as

$$P_{rim} \approx \frac{a}{2} \text{Re} \left\{ \iint E_z(\rho=a) H_x^*(\rho=a) \sin \phi \, d\phi \, dz \right\} - \frac{a}{2} \text{Re} \left\{ \iint E_z(\rho=a) H_y^*(\rho=a) \cos \phi \, d\phi \, dz \right\}. \quad (19)$$

Now, the $E_z$ field is split into three contributions $E_{xz}$, $E_{zy}$, and $E_{zz}$, corresponding respectively to the x-directed surface current, y-directed surface current and the z-directed excitation current. In the same way, we define $H_{xx}$, $H_{xy}$, $H_{xz}$, $H_{yx}$, $H_{yy}$, and $H_{yz}$. The spatial field evaluation at the antenna radius, is carried out by computing with a complex-contour deformation, the relevant Sommerfeld integrals with the appropriate Green's functions. Those integrands converge relatively fast, given the narrow bandwidth of the FBBF spectrum. However, when evaluating the feeding contribution, at (or close to) the feeding layer, the integrand does not converge well. In this case, one should extract an asymptotic term corresponding to the homogeneous medium Green's function while adding this term explicitly in spatial domain. This procedure is explained in the supplementary material.

To efficiently compute azimuthally the rim contributions, the key issue consists of expressing the azimuthal variation of the fields in (19) through Fourier harmonics. Then, using Parseval's identity, one can compute the power in (19) as a single summation over the harmonics. Integration along $z$ is carried out numerically, with a few sampling points (about 4). Details regarding the fields expansion procedure into Fourier Harmonics can be found in the supplementary material as well as in Ref.³¹, Section 7.3.2.

**Ohmic power losses.** This section is devoted to the Ohmic power losses evaluation in the MTS substrate. It is assumed that the losses in the metallization are negligible compared to the substrate loss, which is a good assumption at frequencies lower than 30 GHz²⁹. Considering a substrate with electric conductivity $\sigma$, the Ohmic power loss is given by the third term in (10), where the integration is carried out over the volume just beneath the MTS of radius $a$, over a height corresponding to the substrate thickness $h$. In order to efficiently compute the electric field in the substrate, we have defined a substrate interaction matrix $[Z_G^S(z)]$ between the surface current $J$ and the scattered electric field in a given planar layer in the substrate (fixed height $z$) of radius $a$. The interaction matrix $[Z_G^S(z)]$, as well as the excitation fields vector $[V^E(z)]$, can be computed in terms of FBBFs using a procedure similar to what has been described in²⁷. Therefore, the total electric field at a given height $z$ in the substrate, tested with the complex conjugate of the FBBFs $[E(z)]$ can be calculated as

$$[E(z)] = [Z_G^S(z)] [I] - [V^E(z)]. \quad (20)$$

Based on (20), the coefficients $e_{mn}$ of the total electric field in FBBFs are derived from the tested fields through renormalization. After expanding each component of the electric field into FBBFs and using the orthogonality relation of the FBBFs, the Ohmic power loss in (10) can be rewritten as:
$$P_{\Omega} = \frac{\sigma}{2} \sum_{m,n} K(m, n) \int_{z=-h}^{0} (e_{mn}^{x*} e_{mn}^{x} + e_{mn}^{y*} e_{mn}^{y} + e_{mn}^{z*} e_{mn}^{z}) dz, \quad (21)$$

which means that the Ohmic power losses can be computed extremely fast assuming that the fields expansion coefficients are known. The computation time for the fields coefficients is dominated by the computation time of the underlying substrate matrix $[Z_G^S(z)]$, which requires in practice a few seconds (per sampled layer in the substrate) on a traditional laptop computer.

**Efficiencies under consideration.** From the previous sections, the total power delivered by the feed can be calculated in two different ways i.e

$$P_{tot} = P_{rim} + P_{MTS} + P_{\Omega} = P_{vis} + P_{diff}^{sw} + P_{\Omega}. \quad (22)$$

The first method is based on FBBFs and the second one uses GRBFs. It is worth mentioning that although the second method is implemented with GRBFs, this approach can also be used with any kind of basis functions. Now, we define the following efficiencies.

**Conversion efficiency.** The conversion efficiency $\epsilon_{conv}$ is defined as the ratio between the radiated power $P_{vis}$ and the total power delivered by the feed

$$\epsilon_{conv} = \frac{P_{vis}}{P_{tot}}. \quad (23)$$

Note that the radiated power is computed taking also into account the space waves directly radiated by the feed. Therefore, the conversion efficiency is only affected by the losses and the reflection/diffraction effects at the substrate rim. In contrast to what is defined in$^{29}$, the conversion efficiency is defined on the basis of the total radiated power, i.e. the feed efficiency is not separately treated here.

**Tapering efficiency.** The tapering (aperture) efficiency $\epsilon_{tap}$ is the ratio between the antenna directivity and the maximum directivity that can be achieved considering the antenna area. The maximum directivity is obtained with a uniform illumination of the aperture. For a broadside beam, the tapering efficiency is

$$\epsilon_{tap} = \frac{\lambda^2 D_{MAX}}{4\pi S_1}, \quad (24)$$

where $S_1 = \pi a^2$ is the physical antenna area, $\lambda$ is the free-space wavelength, and $D_{MAX}$ is the maximum antenna directivity. Note that the directivity is computed assuming an infinitely extended substrate. Hence, contrary to the conversion efficiency, diffraction effects originating from the truncation of the substrate are not taken into account in the tapering efficiency.

**Loss factor.** The loss factor $L_{loss}$ quantifies the impact of the losses on the conversion efficiency. It is defined as:

$$L_{loss} = \frac{P_{\Omega}}{P_{tot}}. \quad (25)$$

A small loss factor means that the conversion efficiency is mainly limited by the diffraction effects at the substrate rim.

**Diffraction factor.** The diffraction factor $L_{diff}$ evaluates the impact of the substrate truncation effects on the conversion efficiency. More precisely, it quantifies the power lost in SW's propagating beyond the MTS in the assumption of an infinite grounded dielectric substrate.

$$L_{diff} = \frac{P_{diff}^{sw}}{P_{tot}} = \frac{P_{rim} + P_{MTS} - P_{vis}}{P_{tot}}. \quad (26)$$

A small diffraction factor means that the loss factor plays a major role in the value of conversion efficiency.

**Compound efficiency.** There is no universal figure of merit to define the efficiency of an antenna, since the compromise between the previous metrics will depend on the application of interest. In this paper, we define a "compound efficiency" as the product of the conversion efficiency with the tapering efficiency:

$$\epsilon = \epsilon_{tap} \epsilon_{conv}. \quad (27)$$

This choice is roughly justified by the fact that a zero (or very small) value of the tapering/conversion efficiency automatically leads to a zero (or very small) compound efficiency.

**Numerical and experimental validations.** This section provides numerical and experimental validations of the efficiency computation methods presented in the previous sections.
<figure><img src="image_5.png" /><figcaption>Figure 3. Efficiency (in %) of an isotropic MTS antenna radiating a broadside pencil beam.</figcaption></figure>

**Isotropic MTS.** First, we consider a broadside beam MTS antenna operating at 17 GHz, modeled as an impenetrable scalar impedance given by²⁶

$$Z^{+}(\rho, \phi) = jX_{0}[1 + M \sin(2\pi\rho/p - \phi)], \quad (28)$$

where the average reactance is $X_0 = 0.71 \eta_0$, with $\eta_0$ being the free-space impedance. The period $p$ is defined as $p = \lambda/\sqrt{1 + (X_0/\eta_0)^2}$, with $\lambda$ being the free-space wavelength. $M$ is the modulation index (depth). The MTS radius is fixed to $a = 5.65\lambda$ and the sheet impedance lies on the substrate ROGERS 4350B of relative permittivity $\epsilon_r = 3.66$, and thickness $h = 1.524$ mm. The antenna is fed at its center with a vertical elementary dipole placed in the middle of the substrate ($\rho = 0, z = -h/2$).

Figure 3 illustrates the conversion and tapering efficiencies versus modulation index $M$, computed with the FBBF and the GRBF based methods. The conversion efficiency obtained with the commercial software FEKO³² is reported in the same figure. An excellent agreement is observed between the three approaches, despite being based on completely different basis functions (FEKO uses Rao-Wilton-Glisson (RWG)³³ basis functions) and efficiency computation methods. At $M = 0$, there is a mismatch with FEKO probably due to a different feed modeling. One should note that the FBBF and the GRBF based methods require respectively 2.5 min and 1 min of total computation time on a conventional laptop computer, while FEKO needs 2.4 h. As expected, for relatively low modulation depths, the conversion efficiency is small. However, the conversion efficiency is not zero at $M = 0$ due to the space waves directly radiated by the feed. Moreover, the tapering efficiency is also relatively small for small $M$ because the tapering efficiency takes into account the space waves from the feeder. Indeed, for small modulation depths, the radiation contribution is dominated by those space waves, which are not present at broadside because the feeder is a TM vertical dipole. Increasing the modulation depth, the conversion efficiency improves, but the tapering efficiency reaches a maximum at approximately $M = 0.3$. For higher values of $M$, the tapering efficiency starts to decrease because the uniformity of the tapering illumination rapidly degrades.

**Anisotropic MTS with uniform modulation depth.** A broadside beam anisotropic MTS antenna implementing a sheet impedance modulation is now considered.

$$ \begin{aligned} Z^{\rho\rho}(\rho, \phi) &= jX_0[1 + M(\rho) \cos(\beta_{sw}\rho - \phi)] \\ Z^{\rho\phi}(\rho, \phi) &= jX_0 M(\rho) \sin(\beta_{sw}\rho - \phi) \\ Z^{\phi\phi}(\rho, \phi) &= jX_0[1 - M(\rho) \cos(\beta_{sw}\rho - \phi)], \end{aligned} \quad (29) $$

where $X_0 = -377 \Omega$, and $\beta_{sw}$ is the SW wavenumber supported by a uniform sheet with reactance $X_0$ laying on a grounded substrate of relative permittivity $\epsilon_r = 3$ and thickness $h = 0.762$ mm. The antenna operates at 29.75 GHz with a radius equal to 13.6 cm. Figure 4a shows the efficiencies obtained with the two methods. A very good agreement between the FBBF and the GRBF based methods is observed. The small difference in the aperture efficiency could be traced back to differences in surface current modeling close to the feed region. GRBFs exclude in the current the portion of the MTS occupied by the feed, while FBBFs consider the current distribution on the whole circular domain. A detailed discussion regarding this aspect can be found in²⁷. Anyway, a trend similar to that of the isotropic MTS is observed. However, the tapering efficiency is better in comparison with the isotropic MTS case owing to a better polarization purity¹¹. The overall (compound) efficiency is maximum approximately at $M = 0.4$ and is approximately equal to 40%. That means, a good compromise should be found between the tapering efficiency and the conversion efficiency. Indeed, a relatively low modulation depth will lead to strong
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Efficiency (in %) of an anisotropic MTS antenna radiating a broadside pencil beam. (a) In absence of Ohmic losses. (b) With Ohmic losses in the substrate and using the FBBF based method.</figcaption></figure>

reflection/diffraction effects at the substrate rim which are not predicted by the MoM simulation tools (based on infinite substrate assumption). This diffraction will increase the measured sidelobe level. On the contrary, a high modulation depth will considerably degrade the tapering efficiency of the antenna. Therefore, the methods developed in this paper are useful for the antenna designer to rapidly and accurately optimize the efficiency of the antenna without the need of huge computational resources. Those tools may also be used with MTS antennas designed with the traditional holographic technique<sup>7</sup>; in that case, it provides a rapid way to rigorously check the power balance obtained with holographic theory.

Now, losses are introduced in the substrate through a dissipation factor tan δ = 0.001. The trend of the conversion efficiency, tapering efficiency, diffraction factor and losses factor is reported in Fig. 4b. One can conclude that the losses are very small (about (2 – 3)% of the power delivered by the feeder) and the impact of the losses on the compound efficiency is negligible. This is one of the main advantages of this class of antennas.

*Anisotropic MTS with non uniform modulation depth.* The previous example shows that one can reach 40% compound efficiency by optimizing the value of the uniform modulation depth. However, better performance can be achieved by impressing a non-uniform modulation depth¹⁸. This means that the modulation index increases with respect to the radial coordinate so as to maintain a quasi-uniform aperture field illumination (see Fig. 5a). Such antenna has been designed in¹⁸ using the FO formalism³⁰. Details regarding the design and measurements can be found in¹⁸. The analysis of this MTS provides (with the two methods proposed in this paper) a conversion efficiency of about 95%, which means that almost all the power delivered by the feed is radiated by the antenna. Diffraction at the antenna rim can therefore be neglected. As a result, the directivity analysis with the infinite substrate assumption should provide a good estimate of the measured directivity. Figure 5b compares a cut of the copolar directivity computed using the MoM codes based on FBBF and GRBF with the measurements. The 3D copolar radiation pattern is represented in Fig. 5c–e. One can observe a very good agreement between measurements and simulations. The slight difference in the far sidelobes can be explained by the residual SW power diffracted at the antenna rim, the IBC implementation process, and the modeling of the feeder. The latter is modeled as a simple infinitesimal dipole. Nevertheless, the maximum directivity (about 37.11 dBi) is well predicted by the MoM tools based on GRBF and FBBF, due to the very high conversion efficiency. Lower conversion efficiencies will lead to a higher (in comparison with measurements) predicted maximum directivity. The measured and simulated directivity of the antenna corresponds to about 70% compound efficiency. This efficiency can exceed 80% by minimizing the space waves power directly radiated by the feed. This calls for an efficient SW launcher.

Discussion

Two algorithms have been proposed for the fast evaluation of the efficiency of MTS antennas. Those tools may be used with antenna designs that have been obtained using any type of MTS synthesis, including the traditional holographic method<sup>7</sup>, or the more accurate Flat optics<sup>19</sup>, or the direct inversion method<sup>21</sup>. Whenever the antenna is synthesized, the two formulations presented here provide a rapid algorithm to rigorously check the power balance provided by any of the above mentioned methods. The first formulation is based on the integration of the Poynting vector along a surface enclosing the feed while the second approach relies on the residue evaluation of the electric field at the substrate SW pole. Those algorithms can cope with arbitrary anisotropic surface impedance modulation and provide accurate results in a few minutes. Although based on completely different formalisms and basis functions, those methods provide approximately the same efficiency results and those results have been confirmed with commercial software and experimental data. The method based on the
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Anisotropic MTS antenna with non-uniform modulation index. (a) $X_{\rho\rho}$ reactance at $\phi = 0$. (b) Directivity at $\phi = 0$. The inset disk represents the absolute value of the current distribution on the MTS in log scale. (c) Measured directivity (dBi) in the uv plane. (d) Simulated directivity (dBi) in the uv plane with FBBF. (e) Simulated directivity (dBi) in the uv plane with GRBF. The conversion efficiency of the antenna is about 95% with a 70% compound efficiency.</figcaption></figure>

residue theorem is more analytical and, hence, faster than the Poynting vector based method. Conversely, the Poynting vector based method may be applied to a more general class of problems. Progressive improvement of the antenna compound efficiency is illustrated from isotropic MTS with uniform modulation index (20% efficiency) to anisotropic MTS with optimized (non-uniform) modulation index profile (70% efficiency). To reach and exceed 80% efficiency, one needs to properly design the feeder so as to minimize the space waves directly radiated by the feeder. The study also demonstrated that the Ohmic power losses in modulated MTS antennas are very small (< 10% of the total power) at frequencies lower than 30 GHz. At higher frequencies, one should consider the losses in the metallization. Finally, as in the present paper, the feeder has been modeled as an infinitesimal monopole, the reflection coefficient of the antenna cannot be computed. An accurate estimate of the reflection coefficient should require a full-wave analysis of the MTS including details regarding the feeder.

## Methods

**FBBF and GRBF analysis.** Each FBBF and GRBF is defined according to two parameters *m* and *n* with *m* = 1 ... *M* and *n* = −*N* ... *N* (see equation (2)). The total (x-directed and y-directed) number of basis functions is therefore $2M(2N + 1)$. For the first example, isotropic MTS, we used for FBBF and GRBF $N = 8$ and $M = 46$, which corresponds to 1564 basis functions. In the anisotropic MTS examples, the MTS has been analyzed with FBBF using $N = 10$ and $M = 80$ i.e. 3360 FBBFs, and with GRBF using $N = 10$ and $M = 90$ i.e. 3780 GRBFs.

**FEKO analysis.** FEKO implements the EFIE and allows one to use layered medium Green's functions (FEKO's GF card) as with FBBF and GRBFs. The FEKO simulations have been performed using the impedance boundary condition (IBC) available in this software (SK card and user defined surface impedance). Hence, all the results have been obtained using the same IBC and integral equation, for a fair comparison. Figure 6 shows
<figure><img src="image_8.png" /><figcaption>**Figure 6.** A FEKO simulation setup, each spiral strip of different color stands for a different value of constant surface reactance.</figcaption></figure>

a top view of the simulated structure. Each color strip with spiral shape has a width equal to one tenth of the free space wavelength, and it corresponds to a region with a constant isotropic reactance. Thus, the spiral strips (see Fig. 6) are not used as discrete elements, but to separate the antenna aperture in regions with a constant IBC, i.e., each spiral strip corresponds to a value of surface reactance. Moreover, the structure has been excited using an electric Hertzian dipole (A5 card in FEKO) placed at the origin, in the middle of the dielectric slab. The EFIE currents and radiated fields have been obtained using the standard FEKO solver with 54907 RWG basis functions.

Received: 16 July 2020; Accepted: 5 October 2020

Published online: 15 October 2020

## References

1. Yu, N. et al. Flat optics: controlling wavefronts with optical antenna metasurfaces. *IEEE J. Sel. Top. Quantum Electron.* **19**, 4700423 (2013).

2. Zhang, K. et al. High-efficiency metalenses with switchable functionalities in microwave region. *ACS Appl. Mater. Interfaces* **11**, 28423–28430 (2019).

3. Yueyi, Y. et al. A fully phase-modulated metasurface as an energy-controllable circular polarization router. *Adv. Sci.* **7**, 2001437 (2020).

4. Yueyi, Y. et al. Independent phase modulation for quadruplex polarization channels enabled by chirality-assisted geometric-phase metasurfaces. *Nat. Commun.* **11**, 1–9 (2020).

5. Quarfoth, R. & Sievenpiper, D. Artificial tensor impedance surface waveguides. *IEEE Trans. Antennas Propag.* **61**, 3597–3606 (2013).

6. Mencagli, M., Martini, E., González-Ovejero, D. & Maci, S. Metasurfing by transformation electromagnetics. *IEEE Antennas Wirel. Propag. Lett.* **13**, 1767–1770 (2014).

7. Fong, B. H., Colbrun, J. S., Ottusch, J. J., Visher, J. L. & Sievenpiper, D. F. Scalar and tensor holographic artificial impedance surfaces. *IEEE Trans. Antennas Propag.* **58**, 3212–3221 (2010).

8. Patel, A. M. & Grbic, A. A printed leaky-wave antenna based on a sinusoidally-modulated reactance surface. *IEEE Trans. Antennas Propag.* **59**, 2087–2096 (2011).

9. Maci, S., Minatti, G., Casaletti, M. & Bosiljevac, M. Metasurfing: addressing waves on impenetrable metasurfaces. *IEEE Antennas Wirel. Propag. Lett.* **10**, 1499–1502 (2011).

10. Minatti, G., Maci, S., De Vita, P., Freni, A. & Sabbadini, M. A circularly-polarized isoflux antenna based on anisotropic metasurface. *IEEE Trans. Antennas Propag.* **60**, 4998–5009 (2012).

11. Minatti, G. et al. Modulated metasurface antennas for space: synthesis, analysis and realizations. *IEEE Trans. Antennas Propag.* **63**, 1288–1300 (2014).

12. Li, M., Xiao, S. & Sievenpiper, D. F. Polarization-insensitive holographic surfaces with broadside radiation. *IEEE Trans. Antennas Propag.* **64**, 5272–5280 (2016).

13. Casaletti, M., Smierzchalski, M., Ettorre, M., Sauleau, R. & Capet, N. Polarized beams using scalar metasurfaces. *IEEE Trans. Antennas Propag.* **64**, 3391–3400 (2016).

14. Pandi, S., Balanis, C. A. & Birtcher, C. R. Design of scalar impedance holographic metasurfaces for antenna beam formation with desired polarization. *IEEE Trans. Antennas Propag.* **63**, 3016–3024 (2015).

15. Minatti, G. et al. Metasurface antennas, in Aperture Antennas for Millimeter and Sub-Millimeter Wave Applications (Eds. Springer, 2017).

16. Minatti, G., González-Ovejero, D., Martini, E. & Maci, S. Modulated Metasurface Antennas, in Surface Electromagnetics: With Applications in Antenna, Microwave, and Optical Engineering (Cambridge University Press, Cambridge, 2019).
17. Oliner, A. & Hessel, A. Guided waves on sinusoidally-modulated reactance surfaces. *IRE Trans. Antennas Propag.* **7**, 201–208 (1959).

18. Faenzi, M. et al. Metasurface antennas: new models, applications and realizations. *Sci. Rep.* **9**, 1–14 (2019).

19. Minatti, G., Caminita, F., Martini, E., Sabbadini, M. & Maci, S. Synthesis of modulated-metasurface antennas with amplitude, phase, and polarization control. *IEEE Trans. Antennas Propag.* **64**, 3907–3919 (2016).

20. González-Ovejero, D., Minatti, G., Chattopadhyay, G. & Maci, S. Multibeam by metasurface antennas. *IEEE Trans. Antennas Propag.* **65**, 2923–2930 (2017).

21. Bodehou, M., Craeye, C., Martini, E. & Huynen, I. A quasi-direct method for the surface impedance design of modulated metasurface antennas. *IEEE Trans. Antennas Propag.* **67**, 24–36 (2019).

22. Bodehou, M., Martini, E., Maci, S., Huynen, I. & Craeye, C. Multibeam and beam scanning with modulated metasurfaces. *IEEE Trans. Antennas Propag.* **68**, 1273–1281 (2020).

23. Bodehou, M., Craeye, C. & Huynen, I. Multifrequency band synthesis of modulated metasurface antennas. *IEEE Antennas Wireless Propag. Lett.* **19**, 134–138 (2020).

24. Patel, A. M. & Grbic, A. Modeling and analysis of printed-circuit tensor impedance surfaces. *IEEE Trans. Antennas Propag.* **61**, 211–220 (2013).

25. Francavilla, M. A., Martini, E., Maci, S. & Vecchi, G. On the numerical simulation of metasurfaces with impedance boundary condition integral equations. *IEEE Trans. Antennas Propag.* **63**, 2153–2161 (2015).

26. González-Ovejero, D. & Maci, S. Gaussian ring basis functions for the analysis of modulated metasurface antennas. *IEEE Trans. Antennas Propag.* **63**, 3982–3993 (2015).

27. Bodehou, M., González-Ovejero, D., Craeye, C. & Huynen, I. Method of moments simulation of modulated metasurface antennas with a set of orthogonal entire-domain basis functions. *IEEE Trans. Antennas Propag.* **67**, 1119–1130 (2019).

28. Vegni, L., Cicchetti, R. & Capece, P. Spectral dyadic green's function formulation for planar integrated structures. *IEEE Trans. Antennas Propag.* **36**, 1057–1065 (1998).

29. Minatti, G., Martini, E. & Maci, S. Efficiency of metasurface antennas. *IEEE Trans. Antennas Propag.* **65**, 1532–1541 (2017).

30. Minatti, G., Caminita, F., Martini, E. & Maci, S. Flat optics for leaky-waves on modulated metasurfaces: adiabatic floquet-wave analysis. *IEEE Trans. Antennas Propag.* **64**, 3896–3906 (2016).

31. Bodehou, M. *Modulated metasurface antennas: analysis and direct numerical synthesis*. Ph.D. thesis, UCLouvain (2020).

32. Software, E. & Ltd., S.-S. P. Feko: Field computations involving bodies of arbitrary shape. Online. Available: http://www.feko.info. (2014).

33. Rao, S. M., Wilton, D. R. & Glisson, A. W. Electromagnetic scattering by surfaces of arbitrary shape. *IEEE Trans. Antennas Propag.* **3**, 409–418 (1982).

## Acknowledgements

The authors thank Wave Up Srl (Siena, Italy) for providing measurements of the antenna prototype. This work has been partially funded by Région Wallonne through the LUMINET project under contract number 1810114.

## Author contributions

M.B, C.C and I.H contributed to the development of the Poynting vector based method. D.G, S.M and E.M contributed to the development of the residue based method. M.B implemented the Poynting vector based method. E.M implemented the residue based method. All the authors contributed to the writing of the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-020-74674-w.

Correspondence and requests for materials should be addressed to M.B.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020